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ABSTRACT 

We  show  that  the  total  combinatorial  complexity  of  all  non-convex  cells  in  an 
arrangement  of  n  (possibly  intersecting)  triangles  in  3-space  is 
0{n?'^a.{ny-'^\og^'^n)  and  that  this  bound  is  almost  tight  in  the  worst  case. 
Our  bound  significantly  improves  a  previous  nearly  cubic  bound  of  Pach  and 
Sharir.  We  also  present  a  (nearly)  worst-case  optimal  randomized  algorithm 
for  calculating  a  single  cell  of  the  arrangement  and  an  edternative  less  effi- 
cient, but  still  subcubic  algorithm  for  calculating  all  non-convex  cells,  analyze 
some  special  cases  of  the  problem  where  improved  bounds  (and  faster  algo- 
rithms) can  be  obtained,  and  describe  applications  of  our  results  to  transla- 
tional  motion  planning  for  polyhedra  in  3-space. 

1.   Introduction 

1.1.   Terminology 

Consider  a  collection  G  of  n  (possibly  intersecting)  closed  flat  triangular  (or, 
more  generally,  arbitrary  convex)  objects  Ai,  •  •  ■  ,A„  in  three-dimensional  Euclidean 
space  R^  (see  Figure  1.1).  Let  A  =  A{G)  denote  the  arrangement  of  these  triangles, 
i.e.,  the  subdivision  of  R^  induced  by  them;  thus  A  is  a  decomposition  of  3-space  into 
pairwise  disjoint  connected  cells  of  0,  1,  2,  or  3  dimensions,  where 

(i)     a  3-dimensional  cell  is  a  connected  component  of  R-'  —  UA,-, 

i 

(ii)    a  2-dimensionaI  cell  {di  face)  is  a  connected  component  of  int(Ai)—  U  A ,-  for 

some  i  (where  int(a)  denotes  the  relative  interior  of  A), 

(iii)  a  1-dimensional  cell  (an  edge)  is  either  a  maximal  connected  portion  of  (the  rela- 
tive interior  of)  an  edge  of  some  triangle  A,-,  which  does  not  intersect  any  other 
triangle,  or  a  maximal  connected  portion  of  (the  relative  interior  of)  the  inter- 
section segment  A,- D  A^  of  two  triangles,  which  does  not  meet  a  third  triangle, 
and 
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Figure  1.1.  An  arrangement  of  triangles  in  R  . 

(iv)  a  0-dimensional  cell  (a  vertex)  is  a  vertex  of  a  triangle,  an  intersection  of  an 
edge  of  one  triangle  with  another  triangle,  or  an  intersection  point  of  three  dis- 
tinct triangles. 

To  simplify  the  analysis,  we  henceforth  assume  that  our  triangles  are  placed  in 
general  position  in  space,  meaning  that  no  four  share  a  point,  no  three  intersect  in 
more  than  a  point,  no  two  in  more  than  a  segment,  no  two  edges  of  different  trian- 
gles intersect,  no  triangle  vertex  lies  in  another  triangle,  and  no  triangle  edge  inter- 
sects another  triangle  in  more  than  a  point.  These  assumptions  guarantee  in  particu- 
lar that  each  of  the  above  types  of  cells  have  the  correct  dimensionality.  Since  our 
goal  is  to  analyze  the  combinatorial  complexity  of  (certain  portions  of)  A{G),  this 
assumption  involves  no  real  loss  of  generality,  because  any  degenerate  layout  of  tri- 
angles can  always  be  viewed  as  a  limiting  case  of  layouts  in  general  position,  whose 
combinatorial  complexity  dominates  that  of  the  limit  configuration. 

We  will  use  the  following  additional  terminology.  A  vertex  of  a  triangle  will  be 
called  a  corner;  an  edge  of  a  triangle  will  be  called  an  exposed  edge,  and  the  1- 
dimensional  cells  of  A{G)  into  which  it  is  decomposed  will  be  called  exposed  seg- 
ments. Finally,  the  unqualified  term  cell  will  be  used  to  denote  a  3-dimensional  cell 
of  A(G). 

At  various  points  in  our  analysis  we  will  have  to  consider  planar  arrangements 
of  segments  which  are,  similar  to  the  3-dimensional  case,  partitions  of  a  plane  into 
faces,  edges,  and  vertices  induced  by  a  collection  of  segments  in  that  plane.  As  a 
bounded  face  of  such  an  arrangement  need  not  be  simply  connected,  its  boundary  will 
in  general  consist  of  a  single  outer  component  that  encloses  the  entire  face  and  zero  or 
more  islands  each  enclosing  a  portion  of  the  plane  that  does  not  belong  to  the  face. 
An  unbounded  face  has  a  similar  structure,  except  its  boundary  does  not  have  an 
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outer  component. 

1.2.   Motivation  and  main  results 

Our  goal  is  to  obtain  sharp  upper  and  lower  bounds  on  the  maximum  (worst- 
case)  combinatorial  complexity  of  any  single  cell  C  of  A(G),  that  is  the  number  of 
vertices,  edges,  and  faces  of  A(G)  lying  along  the  boundary  of  C.  Consider  first  the 
following  classification.  A  cell  C  of  A(,G)  is  interesting  if  its  closure  C_contains  at 
least  one  exposed  segment  (it  follows  easily  from  the  definition  that,  if  C  intersects 
an  exposed  segment,  it  fully  contains  it).  All  other  cells  of  A(G)  are  called  dull.  It  is 
easily  checked  that  a  dull  cell  of  A(G)  is  (the  interior  of)  a  convex  polyhedron^ 
bounded  by  (some  of)  the  planes  containing  the  triangles;  its  combinatorial  complex- 
ity is  thus  only  0(n).  Interesting  cells,  on  the  other  hand,  may  have  highly  irregular 
shape;  in  fact,  if  n  is  sufficiently  large,  an  interesting  cell  may  approximate  any  open 
connected  semi-algebraic  set  in  R-',  its  boundary  may  have  arbitrary  genus,  etc.  Let 
C(rt)  be  the  maximum  combinatorial  complexity  of  a  single  (interesting)  cell,  over  all 
possible  arrangements  of  n  triangles  in  general  position.  Simple  examples  show  that 
C(rt)  is  n(n^)  (see  Figure  1.2).  A  more  sophisticated  example  demonstrates  that  C(n) 
is  Xl(/i^a(n))  [PS],  where  a(n)  is  the  extremely  slowly  growing  inverse  Ackermann 
function.  In  fact,  the  combinatorial  complexity  of  the  upper  envelope  of  the  triangles 
in  G  (which  is  a  portion  of  the  boundary  of  the  unbounded  cell  of  A(G),  certainly  an 


Figure  1.2. 


Note  that,  although  in  the  simple  case  of  triangles  (or  arbitrary /7a/  objects)  the  notions  of  dull  and 
interesting  coincide  with  those  of  convex  and  non-convex,  respectively,  this  is  not  necessarily  the  case  for 
more  general  surface  patches. 
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interesting  cell),  is  always  0(n'^a(n))  [PS].  Following  the  results  of  [PS],  as  well  as 
the  work  of  [PSS]  and  [EGS2],  mentioned  below,  on  the  two-dimensional  analog  of 
our  problem,  it  is  natural  to  conjecture  that 

i(n)  =  0(n'^a(n)).  (Conjecture) 

We  have  been  able  to  prove  this  conjecture  only  for  a  few  special  cases,  and  in  gen- 
eral it  is  still  an  open  problem.    The  only  previously  known  non-trivial  upper  bound 

3 — — 
on  5(n),  given  in  [PS],  is  0(n       ^^  );  it  is  proved  using  a  fairly  intricate  analysis 

based  on  extremal  hyper-graph  theory.    (It  is  easy,  by  the  way,  to  show  that  the 

worst-case  combinatorial  complexity  of  the  full  arrangement  /4(G)  is  0(n^).)    Our 

results  (stated  below)  improve  and  generalize  this  bound,  and  we  believe  that  our 

techniques  can  be  extended  to  yield  a  bound  close  to  that  conjectured  above. 

The  main  motivation  (also  discussed  in  [PS],  [EGS])  for  studying  this  problem 
comes  from  motion  planning.  Specifically,  let  B  be  an  arbitrary  polyhedral  object 
bounded  by  k  faces,  edges,  and  vertices  and  free  to  translate  in  three-dimensional 
space  while  avoiding  stationary  polyhedral  obstacles  bounded  by  a  total  of  n  faces, 
edges,  and  vertices.  Using  the  standard  technique  of  expanding  obstacles  in  confi- 
guration space,  it  is  easily  seen  ([PS],  [EGS])  that  the  (3-dimensional)  space  FP  of  all 
free  placements  of  B  can  be  represented  as  a  union  of  certain  (3-dimensional)  cells  of 
A(G),  for  the  collection  G  of  0(kn)  triangles  obtained  as  the  Minkowski  (vector) 
differences  between  faces,  edges,  and  vertices  of  the  obstacles  and  vertices,  edges, 
and  faces  of  B,  respectively.  Moreover,  given  an  initial  placement  zq  of  B,  it  suffices 
to  calculate  only  the  connected  component  of  FP  containing  zq,  because  no  other  por- 
tion of  FP  can  be  reached  from  zq  by  a  collision-free  motion.  Thus,  analysis  of  the 
combinatorial  complexity  and  efficient  calculation  of  a  single  cell  of  A(G)  are  major 
components  in  the  design  and  analysis  of  efficient  algorithms  for  this  translational 
motion  planning  problem. 

Before  presenting  our  results,  let  us  first  review  the  analogous  (and  relatively 
simpler)  situation  in  two  dimensions,  namely  the  case  of  an  arrangement  A  of  n  line 
segments  in  the  plane.  It  is  known  that  the  worst-case  combinatorial  complexity  of  a 
single  face  of  A  is  0(na(n)),  which  is  asymptotically  the  same  as  the  worst-case  com- 
plexity of  the  upper  envelope  of  n  segments  [PSS],  [HS],  [WS].  Also,  the  worst-case 
combinatorial  complexity  of  all  interesting  (i.e.,  non-convex)  faces  of  A  is 
0{n'^'^a(nf^\og^'^n)  [AEGS]  (see  also  rEGS2]);  this  bound  is  almost  tight,  as  the 
complexity  in  question  can  be  n(n '*'•')  [EW].  The  above  upper  bound  is  a  special 
case  of  the  results  in  [EGS2],  [AEGS],  which  show  that  the  complexity  of  any  m  dis- 
tinct faces  in  a  planar  arrangement  of  n  segments  is  0{m^'^~^n^'^'^^^  +  na(n)\ogm), 
for  any  8>0  [EGS2],  but  when  m>n^'^,  this  complexity  is  only 
0(m^'^n^'^a(n)^'^\og^'^n)  [AEGS].  (The  bound  derived  in  [AEGS]  is  actually  some- 
what stronger,  but  our  arguments  will  not  require  its  full  generality.)  It  was  further 
shown  in  [AEGS]  that  the  bound  can  in  faa  be  generalized  to 
Odm^'^p^'^  +  n'*'^p^'^)a(n)'^'^\og^'^n  +  na(n)),  for  all  values  of  m  andp,  where  p  is 
the  number  of  pairs  of  intersecting  segments.  (Note  that  it  yields 
0(pain)^'^\og^'^n  +  na(n))  when  m  =  €Kp)  and  p>n,  which  nearly  coincides  with 
the  actual  ©(p  +  n)  complexity  of  the  entire  arrangement.)  Furthermore,  if  the  m 
faces  are  those  intersected  by  another  (fixed)  segment,  their  total  complexity  is  only 
0{na{n)),    as    has    been    recently    shown    in    [EGPPSS].     These    results    for    two- 
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dimensional  arrangements  are  closely  related  to  our  analysis  and  are  extensively  used 
in  this  paper.   Moreover,  they  clearly  point  toward  the  conjecture  made  above. 

As  already  mentioned,  we  have  not  been  able  to  obtain  a  general  near-quadratic 
upper  bound  for  t,(n).  Instead,  we  consider  here  the  problem  of  analyzing  a  larger 
quantity -the  maximum  total  combinatorial  complexity  <t>(n)  of  all  interesting  cells  in 
any  arrangement  of  n  triangles  in  3-space.  Our  main  result  is: 
Theorem     2.5(c):     <()(«)  =  0(/i'^'^a(n)^'^log'"^/j);     this     is     almost    .tight     since 

We  thus  obtain  a  significant  improvement  over  the  bound  of  [PS],  both  because  our 
bound  is  sharper  and  because  it  applies  to  the  total  complexity  of  all  interesting  cells, 
rather  than  that  of  a  single  one.  In  fact,  we  obtain  a  somewhat  more  general  bound 
that  depends  on  the  number  t  of  pairs  and  the  number  k  of  triples  of  intersecting  tri- 
angles in  G.  Specifically,  in  Theorem  2.5(a)  we  show  that,  if  k^t^n,  the  total 
complexity  of  the  interesting  cells  in  an  arrangement  of  n  triangles  is 

(which  reduces  to  the  previous  bound  by  putting  t=0(n^)  and  k=Oin^)).  More- 
over, our  analysis  can  be  extended  to  obtain  similar  bounds  on  <j)(/z)  for  arrange- 
ments of  arbitrary  convex  flat  objects  (called  convex  plates  in  [PS]),  with  the  possible 
addition  of  a  linear  correction  term  to  accommodate  the  complexity  of  the  "exposed" 
boundary  of  these  plates. 

The  results  outlined  above  generalize  the  0(n^'^\ogn)  upper  bound  of  [EGS3] 
on  the  complexity  of  O(n^)  cells  in  an  arrangement  of  n  planes.  In  fact,  our 
Theorem  2.3  (with  k=0{n^),  t=0{n^))  yields  a  generalization  of  the  bound 
Oim^'^n  log  n  +  n^),  given  in  [EGS3],  for  the  complexity  of  any  m  cells  in  a  plane 
arrangement,  to  arrangements  of  triangles.  However,  our  result  is  valid  only  if  all 
non-convex  cells  are  included  in  the  collection  of  desired  cells,  so  in  our  case  m  has  to 
be  "large",  i.e.  at  least  3/i  +2f  (see  below  for  details  and  explanation  of  this  term). 

We  also  consider  the  task  of  actual  calculation  of  a  single  cell  of  A(G)  and 
present  a  randomized  algorithm,  whose  expected  running  time  is  0{^{n)n  )  for  any 
8>0  (with  the  constant  of  proportionality  depending  on  8).  Moreover,  we  demon- 
strate that  the  same  algorithm  is  faster  when  applied  to  some  special  classes  of 
arrangements.  We  also  describe  an  alternative  algorithm  that  computes  all  interesting 
cells  in  an  arrangement  of  n  triangles  in  expected  time  0{n^'^a.(^n)^'^\o%^^'^n)  and 
0{n'"^a(^n)\ogn)  space,  and  present  a  modification  of  this  algorithm  that  computes 
any  subset  of  cells  in  the  arrangement.  For  details,  refer  to  Section  5.2. 

There  are  two  special  cases  where  both  complexity  bounds  and  algorithmic  per- 
formance can  be  significantly  improved.  One  case  involves  triangles  lying  in  planes 
having  only  a  small  number  /of  orientations  (e.g.,  when  each  triangle  is  parallel  to 
one  of  the  three  coordinate  planes).  In  this  case  we  show  that  the  maximum  com- 
plexity of  a  single  cell  as  well  as  that  of  all  interesting  cells  is  &{n^)  for  a  fixed  /. 
More  precisely,  the  latter  complexity  is  Q(Jh^)  for  f=0{n^'^)  and  the  former  is 
n(n^a(/)).  We  also  provide  an  efficient  0((M  +  n^)/logn)-time  deterministic  algo- 
rithm for  the  calculation  of  any  subset  of  cells  in  such  an  arrangement,  given  a  point 
in  each  cell,  where  M  is  the  total  complexity  of  the  cell(s)  being  computed. 
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The  second  restricted  class  of  arrangements  contains  only  two  types  of  objects: 
(1)  arbitrary  horizontal  convex  plates  (parallel  to  the  Ar>'-plane)  and  (2)  vertical  rec- 
tangles (whose  planes  are  parallel  to  the  z  axis  and  whose  top  and  bottom  sides  are 
parallel  to  the  xy-plane).  In  this  case  we  show  that  the  maximum  complexity  of  a  sin- 
gle cell  is  0(/j^a(«)). 

In  addition,  we  argue  that  our  unmodified  general  algorithm  computes  any  sin- 
gle cell  of  an  arrangement  of  n  triangles,  either  with  only  a  constant  number  of  coor- 
dinates or  in  the  second  special  case  above,  in  expected  time  0(n^'*"^)  (for  any  8>0). 

The  technique  used  to  show  the  0(n^'''a(n)^'-'log^'^n)  upper  bound  on  <j)(n)  is 
relatively  simple.  The  proof  is  based  on  the  results  concerning  the  combinatorial 
complexity  of  many  faces  in  an  arrangement  of  segments  in  the  plane,  summarized 
above.  In  addition,  we  adapt  and  extend  some  of  the  technical  tools  developed  in 
[EGS2],  [EGS3]  for  studying  related  problems.  More  specifically,  supposing  that 
bounds  on  the  complexity  of  all  interesting  cells  in  two  subarrangements  have  been 
obtained,  we  seek  to  establish  a  sharp  relationship  between  these  complexities  and 
the  complexity  of  all  interesting  cells  in  the  arrangement  formed  by  the  union  (i.e., 
overlay)  of  these  two  subarrangements.  Such  relationships  (called  "combination  lem- 
mas") have  been  derived  and  exploited  in  [EGS2],  [AEGS]  for  the  case  of  lines  or 
segments  in  the  plane  and  in  [EGS3]  for  the  case  of  planes  in  3-space.  We  obtain  a 
similar  relationship  for  arrangements  of  triangles  (see  Lemma  2.1  below),  using  a 
proof  technique  that  is  somewhat  different  (and,  we  believe,  simpler)  than  those 
employed  in  [EGS2],  [EGS3]. 

To  analyze  the  time  performance  of  our  main  algorithm  we  prove  an  important 
technical  result  ("the  Slicing  Theorem")  which  may  be  of  independent  interest.  It 
states  roughly  that  any  collection  of  cells  in  an  arrangement  of  triangles  can  be  subdi- 
vided into  tetrahedra  without  substantially  increasing  the  total  complexity  of  this  col- 
lection. 

The  paper  is  organized  as  follows:  In  Section  2  we  obtain  our  main  results  on 
the  combinatorial  complexity  of  all  interesting  cells.  Section  3  derives  the  Slicing 
Theorem.  Sertion  4  discusses  the  special  cases  mentioned  above  and  demonstrates 
the  improved  bounds.  Algorithms  for  computing  single  and  multiple  cells  in  arrange- 
ments of  triangles  are  presented  in  Section  5,  and  a  discussion  of  some  applications 
and  open  problems  is  given  in  Section  6.  Appendix  A  contains  some  basic  facts  con- 
cerning the  topology  of  cells  in  a  three-dimensional  arrangement  of  triangles,  which 
are  needed  for  our  analysis.  Appendix  B  exemplifies  the  usefulness  of  our  tech- 
niques by  presenting  a  simple  alternative  proof  of  (a  variant  of)  the  2-dimensional 
combination  lemma  for  polygonal  regions  used  in  [EGS2]. 

2.   The  complexity  of  all  interesting  cells 

Let  G  =  (Ai,  •  •  •  ,  A„}  be  a  collection  of  triangles  in  3-space  as  in  the  introduction. 
We  determine  the  total  complexity  of  the  interesting  cells  of  A{G)  by  analyzing  a 
more  general  problem: 

Consider  collections  G  of  n  triangles  in  R''  with  at  most  t  intersecting  pairs 
and  at  most  k  intersecting  triples,  and  any  set  P  of  m  ^  2r  +  3n  points,  con- 
taining at  least  one  point  on  each  exposed  segment  of  A(G);  we  wish  to  deter- 
mine the  maximum  complexity  C(m,n,t,k)  of  all  cells  of  A{G)  that  contain 
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at  least  one  of  the  given  points  in  their  interior  or  on  their  boundary,  where 
the  complexity  of  a  cell  containing  multiple  points  is  to  be  counted  only  once, 
and  where  the  maximum  is  taken  over  all  such  collections  G  and  sets  of 
points  P.    (Our  points  will  be  chosen  so  that  each  point  lies  in  the  closure  of 
a  unique  cell.) 
Note  first  that  the  number  of  exposed  segments  in  AiG)  is  2t  +  3n  ^  2(2)  +  3n  =  0(rt  ). 
Indeed,  each  endpoint  of  an  exposed  segment  is  either  a  triangle  corner  or  an  endpoint 
of  the  intersection  segment  of  a  pair  of  triangles.    If  P  contains  only  points  lying  on 
exposed  segments,  one  point  on  each  segment,  the  marked  cells  of  A{G)  (i.e.,  the  cells 
whose  closures  contain  those  points)   are  precisely  all  the   interesting  cells.     Hence, 
C(2t  +  3n,n,t,k)  serves  as  an  upper  bound  on  the  complexity  of  all  interesting  cells  of 
A(G),  and  C(2(l)  +  3n,n,{l),{l))  is  an  upper  bound  on  <t)(/i)  (and,  therefore,  also  on 
C(/i)).    In  the  more  general  problem  formulation,  we  consider  all  interesting  cells  plus 
possibly  some  additional  dull  cells. 

We  use  the  following  divide-and-conquer  approach  to  analyzing  C(m,nj,k).  Parti- 
tion G  into  two  subsets  Gi  and  G2  with  |Gi  |  =  [n/2\,  IG2I  =  \n/2].  Refer  to  the  ele- 
ments of  Gi  as  red  triangles,  and  to  those  of  G2  as  blue.  Obtain  recursively  the  "red" 
marked  cells  ^1,  ■■•,  Rr  of  A(G i),  and  the  "blue"  marked  cells  fli,  •  •  •  ,  S^,  of 
A(G2)  (so  that  the  closure  of  each  /?,  and  each  Bj  contains  at  least  one  point  of  P).  For 
each  point  p,-  €f,  let  R^,  Bbi  be  the  red  and  the  blue  cells  marked  by  p,-,  respectively. 
Let  Ei  be  the  connected  component  of  Rn^Bbt  containing  pi,  i.e.,  the  cell  of  A{G) 
marked  by  pc,  we  refer  to  the  cells  £",-  as  purple  cells.  Let  the  total  complexity  of  all  red 
cells  be  p  and  that  of  all  blue  cells  be  p.  Our  goal  is  to  show  that  the  total  complexity 
of  the  purple  cells  is  at  most 

p  +  P  +  G((m2/3jtl/3  +  ^4/5^1/5)jj(^)2/3iQgl/3^  +  ^^(O). 

We  will  refer  to  this  property,  to  be  proved  below,  as  the  Combination  Lemma  for  an 
arrangement  of  triangles.  It  generalizes  a  similar  combination  lemma  obtained  in 
[EGS3]  for  planes  in  R^. 

2.1.   The  Combination  Lemma 

Before  proceeding  with  the  analysis,  let  us  note  that  the  complexity  of  a  cell  C  is 
simply  the  sum  of  the  number  of  vertices,  edges,  and  faces  of  the  boundary  dC  of  C. 
However,  in  order  to  simplify  the  representation  of  cell  boundaries  for  our  algorithms, 
we  will  count  some  of  these  features  with  multiplicities.  Namely,  imagine  replacing 
each  triangle  A  by  a  "puffy  triangle"  A*  which  is  a  thin  body  bounded  by  two  copies  of 
A,  slightly  curved  outward  from  A  and  seamed  together  at  the  (relative)  boundary  of  A. 
If  the  thickness  of  these  puffy  triangles  is  kept  sufficiently  small,  then  our  general  posi- 
tion assumption  on  the  triangles  in  G  guarantees  that  the  resulting  arrangement  of  the 
puffy  triangles  (more  precisely,  the  common  exterior  of  these  thin  bodies)  maintains  the 
combinatorial  and  topological  structure  of  A{G),  except  that  each  face  of  A{G)  now 
appears  as  two  distinct  faces,  one  on  each  side  of  the  corresponding  puffy  triangle  and, 
similarly,  each  (intersection)  edge  of  A{G)  appears  as  four  distinct  edges,  each  (triple- 
intersection)  vertex  of  A{G)  appears  as  eight  distinct  vertices,  etc.  In  the  subsequent 
analysis,  we  will  follow  this  "puffy  model"  and  take  the  resulting  multiplicities  of  boun- 
dary features  into  account. 
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In  Appendix  A  we  also  show  that  the  number  of  edges  and  vertices  on  the  boun- 
dary dC  of  a  cell  C  is  bounded  by  a  linear  function  of  the  number  of  faces  of  dC  plus  a 
certain  correction  term  which  accounts  for  the  fact  that  the  vertices,  edges,  and  faces 
along  dC  do  not  necessarily  constitute  a  triangulation  of  that  boundary  and  that  C  and/or 
some  of  its  faces  need  not  be  simply  connerted.  It  is  demonstrated  as  well  that  the  total 
sum  of  these  correction  terms,  over  all  cells  in  the  arrangement,  is  only  0{n  +t).  Thus 
the  complexity  of  all  interesting  cells  can  be  estimated  by  simply  bounding  the  total 
number  of  faces  on  the  boundaries  of  such  cells,  and  adding  an  0(n+t)  "overhead" 
term,  which  will  be  absorbed  anyway  in  the  bound  that  we  will  obtain.  In  the 
remainder  of  our  discussion  we  will  use  the  term  "complexity  of  a  cell"  to  refer 
exclusively  to  the  number  of  its  faces. 

Let  us  now  return  to  the  partitioning  of  G  into  red  and  blue  subcollections  G\,  Gi, 
as  above.  Consider  a  resulting  purple  cell  £;  it  is  a  connected  component  of  the  inter- 
section of  a  red  cell  R  and  a  blue  cell  B.  Each  face  of  f  is  a  portion  of  a  (red)  face  of  R 
or  of  a  (blue)  face  of  B.  We  will  separately  bound  the  number  of  blue  faces  and  the 
number  of  red  faces  of  E  (and  sum  these  bounds  over  all  purple  cells).  Consider  red 
faces  first.  Notice  that  the  same  face  of  R  can  yield  a  number  of  purple  faces  that  may 
appear  in  E  as  well  as  in  other  purple  cells.  However,  creation  of  one  purple  face  out  of 
each  red  face  is  already  accounted  for  by  the  overall  "red"  complexity  p,  so  let  us  con- 
centrate on  the  number  of  "extra"  faces  cut  out  of  red  faces  — these  appear  either 
because  a  non-convex  portion  of  the  boundary  of  some  blue  cell  B  may  split  a  red  face  / 
into  a  number  of  disconnected  pieces,  or  because  there  may  be  two  points  of  P  in  (the 
closure  of)  R  belonging  to  different  purple  cells  and  each  of  these  cells  may  have  a  por- 
tion of /on  its  boundary. 

We  start  with  the  red  cells  (i.e.,  those  marked  by  P)  in  A(Gi),  and  proceed  to  con- 
struct the  final  purple  cells  incrementally.  Each  step  of  this  procedure  involves  adding  a 
blue  triangle  A  =  A,  to  the  subarrangement  A,_i  of  all  the  red  and  the  first  /  —  1  blue  tri- 
angles, to  obtain  the  next  subarrangement  A,-  (thus  Aq  =  A{G\)  and  A\„i2]  =A(G)).  A 
"currently  purple"  cell  E^'"^^  (i.e.,  a  cell  of  A,_i  marked  by  P)  may  be  modified  by  the 
addition  of  A  in  one  of  several  ways  — it  may  be  trimmed  but  remain  a  single  cell  with 
the  same  topological  structure,  it  may  remain  a  single  cell  but  be  cut  by  A  in  a  way 
which  changes  the  topology  of  its  boundary,  it  may  be  split  into  two  or  more  subcells, 
each  marked  by  a  point  of  P,  or  its  boundary  may  acquire  a  new  connected  component 
(if  A  is  entirely  contained  in  E^'"^^).  (We  need  not  be  concerned  with  the  last  case 
though,  as  in  such  a  situation  the  introduction  of  A  has  no  influence  on  the  number  of 
red  faces  in  purple  cells.)  We  wish  to  estimate  the  number  of  additional  red  faces 
created  on  the  boundary  of  currently  purple  cells  (in  the  arrangement  A,)  by  the  intro- 
duction of  A.  Summing  the  resulting  bounds  over  all  A,,  and  repeating  the  argument 
for  blue  faces  (starting  with  the  blue  arrangement  and  adding  to  it  red  triangles  one  at  a 
time),  we  will  obtain  an  upper  bound  for  the  purple  complexity  in  terms  of  the  red  and 
blue  complexities. 

Consider  the  intersection  of  A  with  the  current  purple  arrangement  A,_i.  It 
appears  as  an  arrangement  G  =  Q(A)  of  segments  in  the  plane  of  A  (all  clipped  to 
within  A),  each  face  of  which  corresponds  to  a  face  of  some  cell  (actually  to  a  pair  of 
faces  — one  on  each  side  of  A)  in  the  updated  arrangement  A,.  For  convenience,  the  rela- 
tive boundary  3A  of  A  is  also  added  to  Q.  Let  r,  be  the  number  of  segments  (other  than 
the  three  edges  of  A)  participating  in  Q  (in  other  words,  the  number  of  triangles  A^ 
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intersecting  A  for  j  <i)-  Let  ki  be  the  number  of  intersecting  pairs  of  segments  of  Q 
(again,  we  do  not  include  intersections  involving  flA).  As  discussed  in  the  introduction, 
a  bounded  face  of  Q  need  not  be  simply  connected;  its  boundary  consists  of  one  exterior 
connected  component  and  of  zero  or  more  interior  components— so  called  islands. 

We  will  refer  to  a  face  of  Q  incident  to  3A  as  a  "boundary  face."  If  we  erase  the 
three  edges  of  A,  the  boundary  faces  fuse  into  the  single  unbounded  cell  of  the  arrange- 
ment formed  solely  by  the  intersections  of  triangles  Ay  (with  ;<:)  with  A,  so  that  their 
total  complexity,  by  [PSS],  is  at  most  0(/,a(r,)).    (See  Figure  2.1.) 

Observe  that  a  "boundary  face"  /  (such  as  face  /i  in  Figure  2.1)  cannot  split  a 
currently  purple  cell  E  of  A,_i  into  two  subcells,  as /does  not  completely  cut  through  E. 
Nevertheless,  it  can  produce  additional  red  faces.  In  fact,  every  segment  e  on  the  boun- 
dary of /that  does  not  touch  dA  represents  a  cut  of  some  purple  face  /'  of  Ai-i  (note 
that  /'  need  not  be  a  red  face,  but  only  an  earlier  added  blue  face;  also,  /'  need  not  be 
cut  in  two  by  e,  but  just  have  e  reduce  the  number  of  its  islands).  In  other  words,  the 
total  number  of  additional  red  faces  produced  by  all  boundary  faces  on  A  is  no  larger 
than  the  total  complexity  of  all  those  boundary  faces,  i.e.    0(r,a(f,)). 

Let  us  now  consider  an  internal  face  /of  Q(A)  (such  as  faces  /2,  /s  in  Figure  2.1). 
Recall  that  we  restrict  our  attention  to  faces  /which  are  contained  in  a  currently  purple 
cell.  Since  /is  internal,  it  must  either  cut  such  a  cell  E  into  two  subcells  Ei,  E2  or 
change  the  topology  of  the  boundary  of  E  without  splitting  E.  In  the  former  case,  if  one 
of  the  two  resulting  subcells  (say  Ei)  does  not  contain  points  of  P,  no  additional  red 
faces  are  produced  because,  by  our  choice  of  P,  Ei  must  be  dull,  i.e.  convex  (for  other- 
wise its  closure  would  have  to  contain  an  exposed  segment,  and  thus  also  a  point  of  P). 
In  particular,  /  itself  is  convex  and  thus  cuts  each  face  of  E  at  most  once.  Hence  each 
face  of  E  yields  at  most  one  face  of  £2.  as  asserted.   If,  on  the  other  hand,  both  £1  and 


y^-"-^  A       / 

I 

Figure  2.1.  The  arrangement  (2(A)  with  three  'purple'  faces. 
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E2  contain  points  of  P,  the  number  of  new  red  faces  created  is  at  most  proportional  to 
the  complexity  of /(by  an  argument  analogous  to  the  case  of  a  boundary  face).  Let  the 
number  of  such  internal  "splitting"  faces  of  Q(A,)  be  5,-.  Finally,  there  is  the  case  when 
introdurtion  of  an  internal  face  /does  not  cut  E  into  two  subcells,  but  only  changes  the 
topology  of  its  boundary.  Again  the  number  of  additional  red  faces  being  created  is  at 
most  the  number  of  edges  of/.  Let  the  number  of  such  "cutting-but-not-splitting"  faces 
of  Q(^i)  be  c,-.  Hence  the  total  number  of  additional  red  faces  created  by  A,  is 
0(/,a(r,))  plus  the  complexity  of  5,-  +  c,  faces  in  a  planar  arrangement  of  r^  segments  of 
which  it,  pairs  actually  intersect.   The  latter  complexity,  by  [AEGS],  is 

Oaisi  +  cd^'h}'^  +  tf'^k}'^)aiu)^'^\og''hi  +  r,a(r,)).  (*) 

rn/21 

Before  proceeding  with  the  analysis,  let  us  observe  that    ^  ^  — ^  since  in  the  incremen- 

1  =  1 

tal  construction  the  segment  A,- DA.,  if  non-empty,  appears  at  most  once,  namely  in 

rn/21 
G(Aniax{ij})-    Similarly,    2^  ki^sk.    Summing  (*)  and  the  contributions  of  the  boundary 

1  =  1 
faces  over  all  blue  triangles  A,-  we  obtain 

o  [^f  ^  \((si + c,)^'hr' + tf"  kr')c.(uf 'Hog'% + tMti)] 


(  =  1 


rn/21 


=  o\'Z(si+cif"icr'  +  ^trk} 


4/5 1,1/5 


1  =  1 


a(t)^'^\og^'^t  +  0(ta(t)). 


which,  by  Holder's  inequality,  is 


O 


rn/21 

1  =  1 


=  o 


2/3 
X 

■rn/21 

1/3 

-1- 

■r'./2i 

1  =  1 

1  =  1 

rn/21  r''/2l 

1=1         1=1 


2/3 


4/5 


.1/3^,4/5^1/5 


r'./2i 

S    ^1- 

i=  1 


l/5\ 


v2/3i 


.1/3 


a(tr'^\og"^t  +  Oita{t)) 


a(r)2'3iog 


1/3, 


+  0(ta(t)). 


Since  P  contains  m  points,  it  is  impossible  to  make  more  than  m  —  \  cuts,  each  of  which 
splits  a  cell  containing  more  than  one  point  of  P  into  two  subcells  each  containing  at 
least  one  point  of  P.  Hence  ^5,  ^  m  - 1.    As  to  ^^i.  we  prove  in  Appendix  A  the  fol- 

1  1 

lowing 

Proposition  A. 2:  In  an  incremental  construction  of  an  arrangement  of  triangles  with  t 
intersecting  pairs,  as  above,  the  total  number  of  faces  that  cut  cells  without  splitting 
them  is  at  most  0(t). 

Recall  that  by  assumption  m  >2r  + 3/1,  so  25,- +  2c,  ^  wi  +  0(r)  =  0(m).    Putting 


everything  together,  and  repeating  the  argument  for  the  blue  faces,  we  see  that  the  total 
increase  in  the  number  of  faces  of  the  purple  cells  is  at  most 
Oam^'^k^'^  +  t'^'^k^'^)a(t)^'Hog^'^t  +  ra(0).  Notice  that,  if  it  >  r,  this  bound  reduces  to 
^(^2/3^i/3jj(,)2/3,Qgi/3,)     We  have  thus  shown: 

Lemma  2.1  (Combination  Lemma  for  Arrangements  of  Triangles):  Let  G  be  a  set  of  n  tri- 
angles in  R-',  let  t  and  k  be  the  number  of  pairs  and  of  triples  of  triangles  of  G  which 
intersect,  and  let  P  be  a  set  of  m  points,  with  at  least  one  point  of  P  lying  on  each 


Triangles  in  space 


-  10  - 


October  2],  1988 


exposed  segment  of  the  arrangement  AiG).  Partition  G  into  two  sets  Gi,  G2,  and 
denote  by  p  (respectively,  P)  the  complexity  (i.e.,  the  number  of  faces)  of  all  cells  of 
A(Gi)  (respectively,  A(G2))  marked  by  P.  Then  the  total  complexity  of  the  cells  of 
A(G)  marked  by  P  is  at  most 

p  +  ^  +  Oam^'^k^'^  +  t'^'^k^'^)oL(tf^\og^'^t  +  tait)). 

The  recurrence  relation  for  C(m,n,t,k)  that  we  want  to  develop  next  also  depends 
on  the  following  graph-theoretic  lemma.  For  a  hypergraph  H,  let  e{H)  be  the  number 
of  edges  in  H. 

Lemma  2.2:  Given  an  arbitrary  3-uniform  hypergraph  H  on  n  vertices,  let  {A,B}  be  a 
partition  of  its  vertices  with  \A  |=ln/2j  and  \B  |=[rt/2l.  Let  H^  (resp.  H2)  be  the 
sub-hypergraph  of  H  spanned  by  A  (resp.  B).   Then  for  some  choice  of  A  and  B 

e(Hi)  +  e(H2)<jeiH). 

Proof:  For  simplicity,  assume  for  the  remainder  of  the  proof  that  n  is  even.  Let  A  be  a 
randomly  selected  set  of  nil  vertices  of  H  (such  that  all  sets  of  size  nil  are  chosen  with 
equal  probability).   Let  B  be  the  set  of  remaining  vertices.   For  an  edge  e  in  H, 

n       r        •         rr     -•  n       r        .         rr      i  ^tlll—'i'        ^        1 

In  particular,  the  expected  number  of  edges  in  Hi  and  H2  together  is 

e\p[e(Hi)  +  e(H2)]  =  e {H)x(Pr  [e  in H i]  +  Pr[e  in H 2])  <  jeiH), 

implying  the  existence  of  a  choice  of  A  and  B  with  eiHi)  +  e(H2)<~e(H),  as  asserted. 

Using  the  above  result  in  conjunction  with  the  Combination  Lemma,  we  deduce: 

Theorem  2.3:  If  G  is  a  set  of  n  triangles  in  R^  with  t  intersecting  pairs  and  k  intersecting 
triples  and  P  is  a  set  of  m  points,  with  at  least  one  point  of  P  lying  on  each  exposed  seg- 
ment of  the  arrangement  A  =  A(G),  the  total  complexity  of  all  cells  of  A  marked  by  P  is 

C(m,n,t,k)  =  0(m^'^k^'^a(tf^\og^'^tlogn+t'^'^k^'^a(t)^'^\og^'^t  +  n  +  ta(t)logn). 

Proof:  Observe  that  the  total  complexity  of  an  arbitrary  arrangement  of  triangles  is  pro- 
portional to  the  number  of  its  vertices,  which  is  easily  seen  to  be  0(n  +  t+k).  Thus 
C(m,n,t,k)  must  satisfy  the  following  recurrence  relation: 

a(n+t+k),  if  m  >  n+t+k, 

C{m,nll,t',k')  +  C(m,nll,t",k") 

+  bam^'^k^'^  +  t'^'^k^^^)a(t)^'^\og^'^t  +  ta(t)),        otherwise. 


C(m,n,t,k):< 


where  a  and  b  are  constants,  and  t' ,  k' ,  t" ,  k"  are  the  numbers  of  pairs  and  triples  of 
intersecting  triangles  within  each  of  the  two  subcollections  G\,  G2,  respectively. 
Lemma  2.2  applied  to  the  triple-intersection  hypergraph  of  the  triangles  (i.e.,  the  hyper- 
graph whose  nodes  are  triangles  and  which  contains  an  edge  {A,-,Aj,A;t}  whenever 
A,nA^nAfc£0  for  distinct  i,  ;,  k)  implies  that  there  is  a  partitioning  of  G  into  two 
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subsets  Gi  and  G2  each  of  size  roughly  n/2  such  that  k'  +k"  ^Vtk.  Also  trivially 
t'+t"<t. 

We  claim  that 

C(m,n,t,k)  <  d(m^'^k^'^\ogn+t^'^k^'^)a{t)^'^]og^'^t  +  e n  +ftait)\ogn,  (1) 

for  some  constants  d,  e,  f,  depending  on  a  and  b.  First,  note  that  if  m^n  +  t  +  k,  (1)  is 
trivially  satisfied  provided  d,e,f>a.  Let  us  turn  our  attention  to  the  general  case.  For 
the  recurrence  to  be  satisfied,  the  following  inequality  must  hold: 

d(m^'^k^'^\ogn  +  t'^'^k'^'^)ait)^'^]og^'^t  +  e  n  +fta(t)\ogn  > 

d(m^'^k'^'Hogin/2)+t"^'h''^'^)a(t'f^\og^'^t'  +  en/2+ft'ait')\ogin/2) 

dim^'^k"'^'^\ogin/2)  +  t'  '^'^k' '  ^'^)a(t"f'Hog^'^t' '  +en/2+  fi"a(t' ')  log  (n/2) 

The  terms  linear  in  n  cancel  and  t',t"^t,  so  it  is  sufficient  to  ensure  that  the  following 
three  inequalities  hold: 

^Jt^'^logn  >  dk''^'^\og(n/2)  +  dk"'^'^\og(n/2)  +  bk^'^ ,  (2) 

dt^'^k^'^  >  dt'^^k''^'^  +  dt"^'^"^'^  +  bt^'^k^'\  (3) 

and 

/ra(r)  logn  >  //'a(r')  log(n/2)  +  ft"oi{t' ')  log(«/2)  +  fem(r) .  (4) 

First  of  all,  by  Holder's  inequality, 

dt'^'^k'^'^  +  dt"^'^k"'^'^  +  bt'^'^k^'^  <  d{t'  +  t")^'\k'  +  k")^'^  +  fer'^'^A:!'^  s 

^^4/5(  U)l/5  +  ^^4/5^1/5  <  (2-2'5j  +  ^)r^'5^1'^ 

4 
Thus  a  sufficiently  large  J  >  fo  ensures  that  (3)  holds.   Now 

ft'a{t')\og{nl2)  +/ir"a(r")log(n/2)  +  Z?ra(r)  < 

yir'a(01og(n/2)  +yir"a(r)log(rt/2)  +  fera(r)  < 

fta{t)\oginl2)  +  fe/a(0  <  /ira(r)  logn  +  {b  -f)ta(t). 
Thus  (4)  holds,  provided  that/>fc. 

Turning  to  (2),  an  application  of  Holder's  inequality  yields 

^,1/3   +  ^,,1/3  <  (A:'+it")^'322/3  <  (1a:)1'322/3  =  ^1/3 

4 

In  particular, 

dk'^'Hogin/2)  +  dk"^'^\og(n/2)  4-  Z^/t^'^  < 

^)t^'3log(«/2)  +  feit^'^  :s  dk^'hogn  +  ib-d)k^'^, 
ensuring  that  (2)  holds  for  any  d>b,  and  thus  completing  the  proof.   □ 
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For  developing  a  lower  bound,  it  will  be  convenient  to  have  the  following  graph- 
theoretical  result: 

Lemma  2.4:  Let  H  be  an  undirected  graph  with  E  edges  and  C  3-cycles.  Then 
C  =  0(E^'^). 

Proof:  Immediate  from  the  following  theorem  of  [E]  that  gives  the  maximum  number  of 
complete  subgraphs  on  r  vertices  contained  in  a  graph  with  E  edges  (in  our  case,  r  =  3). 

Let£:  =  (p  +  /,  0<f<5,  thenC,<(p  +  (^lj) 
D 

Having  proved  Theorem  2.3,  we  obtain  our  main  results  concerning  all  interesting 
cells  by  substituting  m  =  3n  +  2t  and  then  observing  that  t  =  0{n^)  and  k  =  O(n^): 

Theorem  2.5:  (a)  The  complexity  of  all  interesting  cells  in  an  arrangement  of  n  triangles 
with  t  intersecting  pairs  and  k  intersecting  triples  is 

Oiin  +  t)^'^k^'^a(tf^\og^'^t\ogn  +  t'^'^k^'^a(t)^'^\og^'^t  +  n+  ra(r)logn), 

which  reduces  to  0(t'^'^k^''^a(t)^''\og^'^t),  if  k^t^n.  Moreover,  for  any  k>t>:n  such 
that  k  =  0(P'^)  there  exists  an  arrangement  of  n  triangles  with  0(t)  intersecting  pairs 
and  O^k)  intersecting  triples  so  that  the  complexity  of  all  interesting  cells  is 

n(t^'^k^'^  +  ta(t)-\-n). 

(b)  The  said  complexity  is  also  0{(n  +  t)^'^t^'^a(t)^'^\og^'^t\ogn  +n),  regardless  of  the 
value  of  k;  this  is  bounded  by  0(t'^'^ait)^'^\og'^'^t) ,  provided  t>n. 

(c)  In  the  worst  case,  ^(n)  =  0(n^'^a(n)^'^\og'^'-^n).  This  is  almost  tight,  since 
<t)(n)  =  n(/.7'3) 

Proof:  The  upper  bounds  in  (a)  and  (c)  are  immediate  from  the  preceding  analysis, 
while  the  bound  in  (b)  easily  follows  from  that  in  (a)  by  substituting  k  =  0(r  ),  which 
is  a  consequence  of  Lemma  2.4  applied  to  the  intersection  graph  of  the  triangles  (note 
that  a  triple  intersection  of  triangles  is  necessarily  a  3-cycle  in  this  graph,  while  the  con- 
verse is  not  true). 

Turning  to  the  lower  bound  in  (a),  let  us  assume  k^t^n  and  put  a  =  [k/t\  and 
b  =  [t  /k\.  Notice  that  t^  k  hy  assumption,  so  both  a  and  b  are  positive  integers.  To 
obtain  the  lower  bound  in  (a),  consider  a  set  of  a  vertical  rectangles  whose  bottom  edges 
lie  in  the  plane  z  =  0,  whose  top  edges  lie  in  the  plane  z  =  b  +  \,  and  such  that  their  xy 
projections  yield  a  planar  arrangement  of  a  segments  for  which  the  complexity  of  all 
interesting  (i.e.,  non-convex)  faces  is  il(a^''^)  [EW].  Cutting  these  rectangles  by  b  addi- 
tional horizontal  (and  sufficiently  large)  rectangles  lying  in  the  planes  z  =  i,  for 
i=  l,...,b,  is  easily  seen  to  yield  an  arrangement  of  n=a+b  rectangles  (that  can  be 
transformed  into  an  arrangement  of  n  triangles  in  general  position)  in  which  the  total 
complexity  of  all  interesting  cells  is 

Notice  that  the  number  of  pairwise  intersections  in  this  arrangement  is  no  more  than 

7  k^  -       , 

a  +ab  =  —  +  t=0(t)  (as  k  ^t^  by  assumption),  while  the  number  of  triples  of  inter- 
secting triangles  is  at  most  a^b^k.   This  completes  the  construction  of  an  arrangement 
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of  triangles  with  0(t)  intersecting  pairs  and  at  most  k  intersecting  triples  in  which  the 
complexity  of  all  interesting  cells  is  ft{r^'^^^'^).  In  addition,  note  that  the  interesting 
cells  always  have  complexity  Cl(n).  Moreover,  the  construction  in  [PS]  of  an  arrange- 
ment of  n  triangles  in  which  the  unbounded  cell  has  complexity  Cl(n^ain))  can  easily  be 
modified  to  yield  a  similar  construction  with  the  unbounded  cell  having  complexity 
ft(ra(r)).  Putting  the  three  bounds  together,  we  obtain  the  lower  bound  in  (a).  The 
lower  bound  in  (c)  follows  by  substituting  k  =  0(n^)  and  t  =  0(n^).   Q 

Remark:  The  reader  should  note  the  strong  dependency  of  our  bounds  on  the  bounds 
for  planar  arrangements  of  segments.  Any  future  improvements  in  the  two-dimensional 
bounds  (such  as  the  recent  results  in  [AEGS])  will  most  likely  carry  over  to  our  analysis 
and  yield  corresponding  improvements  in  our  bounds.  The  exotic-looking  expression 
a(r)^'''log^'^r  appearing  in  our  bounds  is  a  direct  "carry-over"  from  the  analysis  of 
[AEGS].  Removal  of  this  expression  from  the  bounds  in  [AEGS]  would  also  remove  it 
from  our  bounds.  For  a  more  specific  example,  if  the  complexity  of  m  faces  in  an 
arrangement  of  n  segments  intersecting  in  p  points  should  be  shown  to  be 
0(m^'^p^'''+na(«)),     then    the    bound     in    Theorem     2.5(a)     would     go     down     to 

3.   The  Slicing  Theorem 

In  this  section  we  will  prove  an  auxiliary  result  on  triangulating  cells  in  arrange- 
ments of  triangles.  It  facilitates  a  general  efficient  algorithm  for  computing  a  single  cell 
in  such  an  arrangement  (cf.  Section  4).  A  similar  technique  is  also  used  to  establish 
better  bounds  and  construct  more  efficient  algorithms  for  some  special  classes  of 
arrangements  discussed  below.  In  addition,  our  result  provides  a  partial  answer  to  the 
following  seemingly  simple  question:  given  an  arbitrary  polyhedral  region  AT  in  R^,  can 
it  be  cut  into  a  small  number  of  (pairwise  disjoint)  tetrahedra?  The  two-dimensional 
analog  of  this  problem  has  a  satisfactory  solution  — the  number  of  triangles  needed  to  tri- 
angulate an  arbitrary  planar  polygonal  region  is  proportional  to  its  complexity.  In  three 
dimensions,  the  problem  becomes  much  more  difficult.  Some  negative  results  are  given 
in  [OR,  Chap.  10];  our  result  shows  roughly  that  the  required  number  of  tetrahedra  is 
not  much  larger  than  the  complexity  of  K,  if  AT  is  a  cell  (or  a  collection  of  cells)  in  an 
arrangement  of  triangles.  The  technique  that  we  use  in  the  proof  of  the  Slicing 
Theorem  is  related  to  that  employed  by  Chazelle  and  Palios  [CP]  in  obtaining  similar 
results  for  general  polyhedra. 

Theorem  3.1  (Slicing  Theorem):  Let  AT  be  a  collection  of  cells  in  an  arrangement  of  n  tri- 
angles in  3-space,  with  a  total  of  h  faces.  Then  K  can  be  decomposed  ("triangulated") 
into  O(n^ain)  +  h)  tetrahedra  having  pairwise  disjoint  interiors. 

Proof:  We  construct  the  triangulation  incrementally,  by  adding  new  vertical  faces  (to 
which  we  refer  as  walls)  emanating  from  exposed  edges  of  the  given  triangles.  These 
vertical  walls  will  collectively  decompose  the  cells  of  K  into  convex  subcells.  Once  this 
is  done,  the  final  triangulation  is  easily  obtained  by  triangulating  the  boundary  of  each 
convex  subcell  C  and  connecting  each  resulting  triangle  to  some  fixed  interior  point  of 
C.  The  number  of  final  tetrahedra  is  clearly  proportional  to  the  total  complexity  of  all 
subcells  into  which  the  walls  partition  K,  which  we  now  proceed  to  estimate. 

First  of  all,  we  assume  that  all  cells  of  K  are  interesting;  the  above  observations 
imply  that  this  involves  no  real  loss  of  generality.  Let  e\,  ei,  •  •  •  ,  ^3„  be  the  exposed 
edges  of  the  triangles.    For  each  ei  in  turn  we  add  vertical  walls  emanating  from  it,  as 
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follows.  Suppose  this  has  already  been  done  for  all  ej,  j<i.  Let  p,-  be  the  vertical 
plane  containing  g,-,  and  let  V,-  be  the  collection  of  segments  in  pi,  each  of  which  is  either 
the  intersection  of  p,-  with  some  triangle  or  a  connected  component  of  the  intersection  of 
Pi  with  some  previously  erected  wall.  Then  the  walls  added  at  the  current  (/th)  stage 
are  simply  all  the  horizon  faces  of  e,-  (i.e.,  faces  incident  to  e,)  in  the  planar  arrangement 
A(Vi)  which  are  contained  in  K.  See  Figure  3.1  for  an  illustration.  (For  convenience  of 
exposition,  we  will  need  to  assume  that  no  exposed  edge  is  vertical  and  that  pi^Pj 
whenever  i^^j;  this  can  always  be  enforced  by  an  appropriate  rotation.) 

It  is  easy  to  check  inductively  that  the  addition  of  all  these  walls  results  in  a  convex 
decomposition  of  K.  Indeed,  the  walls  added  at  the  /th  stage  remove  all  non-convex 
edges  occurring  on  «,-,  and  no  new  non-convexities  are  created. 

We  next  analyze  the  total  complexity  (i.e.  the  number  of  faces)  of  the  resulting 
decomposition  of  K.  It  is  easily  seen  that  this  complexity  is  A  +  0(q),  where  q  is  the 
total  complexity  (i.e.,  total  number  of  edges)  of  all  the  vertical  walls.  A  major  obstacle 
in  estimating  q  is  that  we  have  no  a  priori  linear  bound  on  the  number  of  segments  in 
each  V,  — each  triangle  contributes  at  most  one  segment  to  V,-,  but  a  vertical  wall  may 
intersect  p,-  in  many  segments  and  there  are  many  vertical  walls  erected  from  each  ej. 
We  overcome  this  difficulty  as  follows: 

Fix  an  exposed  edge  e,-,  and  let  V^  be  the  collection  of  the  (at  most  n)  segments 
formed  by  intersecting  p,-  with  the  given  triangles.  By  the  result  of  [EGPPSS],  the  total 
complexity  of  all  horizon  faces  of  e,  in  A(V?)  is  0(na(n)).  Let  H  be  the  collection  of 
all  horizon  faces  (within  K)  of  g,-  in  A(V,)  (namely  the  vertical  walls  erected  from  e,). 


Figure  3.1. 

The  arrangement  A(y,)  with  three  previously  erected  sets  of  walls  intersecting  it. 
The  shaded  regions  are  the  new  walls.  Vertical  segments  are  marked  with  their  type. 
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Clearly,  each  face  in  //  is  a  subface  of  some  horizon  face  in  A(Vf).  Note  that,  for  each 
j  <  i,  all  vertical  walls  erected  from  Cj  intersect  p,-  in  a  collection  of  segments  lying  in 
the  vertical  line  Pir\pj  (see  Figure  3.1)  so  that,  in  particular,  at  most  one  of  them  cuts 
g/.  It  follows  that  the  horizon  faces  in  A(V,)  meet  ei  in  at  most  4n  —  \  edges,  because  the 
segments  of  V,  cut  ei  in  at  most  An— 2  points  (there  are  at  most  3/1  —  1  planes  pj  with 
j<i  and  at  most  n  —  1  triangles  not  containing  e,).  Since  each  such  edge  bounds  at  most 
two  faces,  we  conclude  that  there  are  at  most  8n— 2  horizon  faces  in  A(Vi).  We  will 
subdivide  each  segment  of  V,  —  V?  into  subsegments,  cutting  it  at  its  points  of  intersec- 
tion with  segments  of  V?,  and  work  with  the  refined  vertical  segments  from  now  on. 
For  simplicity  assume  that  all  segments  are  finite;  however,  the  argument  that  follows 
can  be  easily  adapted  to  handle  rays  and  full  lines.  Since  each  vertex  of  a  face  in  H  is 
either  a  vertex  of  a  horizon  face  of  e,  in  A(V^)  (of  which  there  are  only  0(na(n))),  or 
an  endpoint  of  some  vertical  segment  in  V,  — Vf,  it  therefore  suffices  to  estimate  the 
total  number  of  such  segments  which  bound  horizon  faces  of  e,-  in  A(V,).  We  now  add 
these  vertical  segments  one  at  a  time,  thus  gradually  transforming  the  horizon  faces  in 
A(Vi*)  into  those  in  A(V,)-  Let  Vj  denote  the  current  collection  of  segments  obtained 
after  adding  j  of  these  segments.  Each  segment  s  in  V,  — Vj  bounding  a  final  horizon 
face  in  A(V,)  is  of  one  of  the  following  types  (note  that,  by  construction,  each  endpoint 
of  a  (finite)  vertical  segment  lies  on  a  segment  of  V^): 

(1)  s  bounds  the  same  horizon  face  in  A{Viu{s})  on  both  of  its  sides, 

(2)  s  lies  on  the  common  boundary  of  two  distinct  horizon  faces  of  A  {Viu{s}),  or 

(3)  s  bounds,  on  one  side,  a  horizon  face  of  A(Viu{s})  and  on  the  other  side  a  face  of 
A(ViU{s})  which  is  no  longer  adjacent  to  «,-. 

Refer  to  Figure  3.1.  Observe  that,  as  a  new  vertical  segment  is  added  to  Vj  to  obtain 
Vi"^^,  the  type  of  the  remaining  segments  may  change;  in  fact,  it  is  easily  verified  that 
the  only  transitions  possible  are  (1)-K2),  (2)-K3),  and  (1)-K3).  We  proceed  to  transform 
the  horizon  faces  of  A(Vf)  into  those  of  A(V,)  by  first  adding  all  segments  of  V,  — Vf 
which  at  the  time  of  addition  have  type  (1).  (As  noted,  no  transition  can  revert  a  seg- 
ment to  type  (1).)  Each  such  segment  reduces  the  number  of  "islands"  (i.e.,  connected 
components  of  the  boundary)  of  some  horizon  face  by  1.  But  the  total  number  of  islands 
in  all  horizon  faces  is  at  most  n,  because  each  segment  of  V?  appears  in  at  most  one 
island,  and  (on  the  assumption  that  no  vertical  segment  is  infinite)  each  island  must  con- 
tain a  segment  of  V^ .  Thus  the  total  number  of  type  (1)  segments  added  is  at  most  n. 
A  similar  argument  shows  that  the  number  of  type  (2)  segments,  which  we  add  after  the 
type  (1)  segments,  is  also  0(n),  because  addition  of  each  such  segment  increases  the 
number  of  horizon  faces  by  1,  and  there  are  only  Oin)  such  faces  in  the  final  A(V,). 
The  remaining  vertical  segments  all  have  type  (3)  and  will  never  change  their  type 
again.  Let  V*  be  the  union  of  Vf  with  all  type  (1)  and  type  (2)  vertical  segments.  Since 
|V*  I  =  0{n),  it  is  still  the  case  that  the  complexity  of  all  horizon  faces  of  e,-  in  A(V*)  is 
0(na(n)).  Now  add  the  vertical  segments  of  type  (3).  Each  such  segment  s,  in  its  turn, 
chops  off  some  portion  f,  of  a  horizon  face,  which  becomes  disconnected  from  «,-. 
Moreover,  since  we  add  only  segments  that  actually  appear  on  the  boundaries  of  the 
final  horizon  faces  in  A(V,)  and  vertical  segments  are  pairwise  disjoint,  each  portion  f, 
is  chopped  off  only  once  (namely,  these  portions  have  pairwise  disjoint  interiors)  and, 
since  f,  is  a  polygonal  region,  it  must  contain  a  vertex  of  a  horizon  face  in  A(V')  other 
than  either  of  the  endpoints  of  s.  Hence  the  number  of  portions  f,,  and  thus  also  the 
number  of  type  (3)  segments  s,  is  Oina(n)).  Repeating  the  analysis  for  each  «,-,  we 
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conclude  that  the  total  number  of  faces  in  the  decomposition  of  K  is  h  +  0(n^a(n)),  as 
asserted.  As  discussed  above,  a  final  triangulation  of  this  decomposition  completes  the 
proof  of  the  theorem.  □ 

Corollary  3.2:  Let  ^  be  a  collection  of  p  cells  in  an  arrangement  of  n  triangles  in  R^, 
having  a  total  of  h  faces.  Then  K  can  be  decomposed  into  p  +  0(n  )  convex  polyhedra 
with  pairwise  disjoint  interiors,  such  that  their  total  complexity  is  only  h  +  0(n^a(n)). 

Proof.  It  is  sufficient  to  demonstrate  that  in  the  above  construction  the  number  of 
resulting  subcells  does  not  exceed  p  by  more  than  O(n^),  since  the  second  claim  follows 
immediately  from  the  proof  of  the  Slicing  Theorem.  Recall  that  a  vertical  wall  w 
erected  at  a  given  step  of  the  construction  is  a  face  in  the  arrangement  induced  in  its 
plane  by  intersections  with  all  triangles  and  all  vertical  walls  erected  in  previous  steps. 
Hence  w  cuts  the  subcell  C  in  which  it  is  contained  into  at  most  two  new  subcells  (note 
that  it  is  possible  for  w  to  change  the  topology  of  the  boundary  of  C  without  cutting  C 
into  two  subcells).  Hence  introduction  of  w  increases  the  number  of  subcells  by  at  most 
1.  However,  there  are  only  0(n)  vertical  walls  erected  from  each  exposed  edge,  thus 
the  total  number  of  convex  polyhedra  obtained  by  the  construction  does  not  exceed  p  by 
more  than  0{n^).  □ 

Corollary  3.3:  The  collection  K  of  all  interesting  cells  in  an  arrangement  of  n  triangles 
in  3-space  can  be  decomposed  into  0{n  )  convex  polyhedra  with  pairwise  disjoint  interi- 
ors, such  that  their  total  complexity  is  only  h  +  0{n^a{n)),  where  h  is  the  complexity  of 
K. 

Proof:  Immediate  from  Corollary  3.2,  since  there  are  at  most  as  many  interesting  cells  as 
there  are  exposed  segments,  i.e.  O(n^).  □ 

Remarks:  (1)  While  the  Slicing  Theorem  is  intuitively  plausible,  the  proof  is  not  at  all 
trivial.  An  open  problem  is  to  generalize  the  theorem  to  the  case  of  an  arbitrary  3- 
dimensional  semi-algebraic  set  K,  namely  to  show  that  K  can  be  decomposed  into  a  col- 
lection of  simple,  connected  cells  (e.g.,  such  as  in  Collins'  cylindrical  algebraic  decom- 
position [Co]),  whose  cardinality  is  roughly  the  same  as  the  combinatorial  complexity  of 
dK.  A  partial  affirmative  result  of  this  sort  is  indicated  in  [CEGSW]  for  the  collection  K 
of  all  3-cells  in  an  arrangement  of  spheres;  the  analysis  given  there  appears  to  be  exten- 
dible to  arbitrary  algebraic  surfaces.  A  related  open  problem  is  to  obtain  generaliza- 
tions of  the  Slicing  Theorem  to  arrangements  of  simplices  in  higher  dimensions. 

(2)  Observe  that  the  key  point  of  the  proof  of  the  Slicing  Theorem  is  the  assertion  that, 
given  a  planar  arrangement  A  of  n  'old'  segments  and  an  unspecified  number  of  'new' 
segments,  the  total  complexity  of  the  faces  of  A  cut  by  a  segment  e  is  still  0{na(^n)), 
provided  the  new  segments  (i)  cross  e  at  most  0{n)  times,  (ii)  do  not  intersect  among 
themselves,  and  (iii)  start  and  end  at  points  of  old  segments  or  "at  infinity".  It  is  easy 
to  verify  that  the  proof  given  above  applies  to  this  slightly  more  general  situation  as 
well. 

(3)  Note  that  even  when  the  number  t  of  pairs  of  intersecting  triangles  in  the  arrange- 
ment is  o(n^),  the  number  of  cells  in  the  resulting  decomposition  of  K  can  still  be  qua- 
dratic (this  is  the  case,  for  example,  in  the  arrangement  of  n  long  and  thin  horizontal 
rectangles,  nil  of  which  lie  in  the  plane  z=0,  while  the  remaining  nil  lie  in  the  plane 
z=l,  as  shown  in  Figure  3.2;  here  in  fact  /=0). 
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Figure  3.2. 

A  collection  of  disjoint  rectangles  in  which  our  "slicing"  produces  O(n^)  polyhedra. 
Solid  rectangles  lie  in  the  plane  z=l,  dashed  ones  in  z  =  0. 


4.   Special  cases 

In  this  section  we  will  discuss  two  restricted  classes  of  arrangements  of  triangles  for 
which  the  worst-case  complexity  of  a  single  cell  (or  of  all  interesting  cells)  is  easier  to 
bound. 

4.1.   Triangles  with  few  orientations 

We  first  consider  arrangements  of  n  triangles  in  which  the  planes  containing  the 
triangles  have  at  most  / distinct  orientations.  (Note  that  in  the  following  analysis /need 
not  be  constant.)  A  (very  restricted)  example  of  such  an  arrangement  would  be  when 
each  triangle  lies  in  a  plane  parallel  to  one  of  the  /=3  coordinate  planes.  In  what  fol- 
lows, we  will  assume  however,  without  loss  of  generality,  that  no  triangle  lies  in  a  plane 
parallel  to  the  z  axis.  Corollary  3.3  yields  a  decomposition  of  all  interesting  cells  in 
such  an  arrangement  into  O(n^)  polyhedra.  Being  convex,  each  of  these  polyhedra  is 
bounded  by  at  most  two  faces  contained  in  planes  having  any  specific  orientation,  so  it  is 
bounded  by  at  most  2/  faces  lying  in  the  original  triangles.  The  remaining  faces  lie  in 
the  vertical  walls  added  in  the  construction  of  the  Slicing  Theorem;  thus  they  do  not  con- 
tribute to  the  complexity  of  the  (undecomposed)  interesting  cells  (except  by  cutting 
other  faces  into  subfaces).  Hence,  the  total  number  of  non-vertical  faces  in  all  polyhe- 
dra lying  in  interesting  cells  is  0(Jh^),  which  is  an  upper  bound  on  the  total  number  of 
faces  in  all  interesting  cells.   This  proves 


Triangles  in  space 


•  18- 


October21.  1988 


Theorem  4.1:  The  complexity  of  all  interesting  cells  in  an  arrangement  of  n  triangles 
lying  in  planes  with  / distinct  orientations  is  at  most  0(Jh^). 

As  for  a  corresponding  lower  bound,  it  is  easy  to  obtain  a  variant  of  the  construc- 
tion of  [EW]  for  a  lower  bound  on  the  complexity  of  many  faces  in  arrangements  of 
lines  which  involves/^  segments  in  the  plane,  having  only  /distinct  orientations,  such 
that  the  complexity  of  all  interesting  (i.e.,  non-convex)  faces  of  their  arrangement  is 
ft(/*).  Assume  f^inll)^'^.  Placing  \nl2f^\  disjoint  translates  of  such  an  arrangement 
side  by  side  in  the  plane  produces  an  arrangement  of  at  most  nil  segments  having /dis- 
tinct orientations,  in  which  the  complexity  of  all  interesting  faces  is 
rt/2/-'  X  n(/*)  =  nc/h).  Now  consider  a  3-dimensional  arrangement  of  nil  vertical  rec- 
tangles whose  lower  edges  lie  in  the  plane  z  =  0,  whose  upper  edges  lie  in  the  plane 
r  =  n/2+l,  and  whose  vertical  projections  induce  the  planar  arrangement  of  nil  seg- 
ments described  above,  together  with  nil  large  horizontal  rectangles  lying  in  the  planes 
z  =  J,  for  z  =  1,  •  •  •  ,nll.  It  is  easy  to  see  that  the  total  complexity  of  the  interesting 
cells  in  this  arrangement  of  n  rectangles  is  at  least  n  x  €l(Jn)  =  iUJn^),  but  only  f+  1  dis- 
tinct plane  orientations  were  used.  Moreover,  the  rectangles  can  be  easily  replaced  by 
triangles  without  reducing  the  combinatorial  complexity  of  the  interesting  cells.  Thus 
we  have  shown 

Theorem  4.2:  The  worst-case  complexity  of  all  interesting  cells  in  an  arrangement  of  n 
triangles  lying  in  planes  with  /distinct  orientations  is  at  least  il(Jh^)  for /=  Oin^'^). 

A  similar  argument,  based  on  the  construction  of  a  planar  arrangement  of  n  seg- 
ments in  which  the  unbounded  face  has  complexity  ft(/ia(rt))  [WS],  gives 

Theorem  4.3:  The  worst-case  complexity  of  a  single  cell  in  an  arrangement  of  n  trian- 
gles lying  in  planes  with  / distinct  orientations  is  at  least  Cl(n  a(f)). 

Corollary  4.4:  In  an  arrangement  of  n  triangles  lying  in  planes  with  a  constant  number 
of  distinct  orientations,  the  worst-case  complexity  of  any  single  cell  (or  of  all  interesting 
cells)  is  e(n^). 

4.2.   The  horizontal-vertical  case 

The  second  special  case  concerns  arrangements  where  only  two  types  of  convex 
plates  are  allowed:  arbitrary  (convex)  horizontal  polygons  (i.e.,  parallel  to  the  ry-plane) 
and  vertical  rectangles  (i.e.,  rectangles  with  two  sides  parallel  to  the  z  axis).  Let  n  be 
the  total  number  of  plates  (both  horizontal  and  vertical)  in  the  arrangement  and  let  q  be 
the  total  complexity  of  all  horizontal  plates.  The  lower  bound  construction  of  Theorem 
2.5  applies  in  this  case,  showing  that  the  total  complexity  of  all  interesting  cells  in  this 
restricted  class  of  arrangements  is  still  fl(n^'^).  We  now  use  a  variant  of  the  decomposi- 
tion given  by  the  Slicing  Theorem  to  prove  a  tight  Q(n^a(n)  +  q)  bound  on  the  worst- 
case  complexity  of  a  single  cell  in  such  an  arrangement.  Recall  that  the  only  type  of 
non-convexity  that  occurs  in  an  interesting  cell  is  along  an  exposed  segment.  We  distin- 
guish two  types  of  interesting  cells:  (1)  those  whose  boundary  contains  only  vertical 
exposed  segments  (i.e.,  those  parallel  to  the  z  axis)  and  (2)  those  whose  boundary  con- 
tains horizontal  exposed  segments  (and  possibly  some  vertical  ones).  We  will  argue  that 
the  complexity  of  a  single  cell  of  type  (1)  is  Oina(n)),  while  the  total  complexity  of 
type  (2)  cells  is  0(n^a(n)  +  q),  concluding  that  the  complexity  of  any  one  cell  in  such  an 
arrangement  is  0(n^a(n)  +  q). 

Consider  a  cell  C  of  type  (1).  Clearly,  it  cannot  be  unbounded.  Moreover,  the 
absence  of  horizontal  exposed  segments  on  dC  implies  that  C  is  in  fact  a  right  vertical 

Triangles  in  space  •  19  -  October  21,  1988 


prism  whose  height  is  a  maximal  vertical  open  segment  extending  from  a  "floor  plate" 
A/  to  a  "ceiling  plate"  A^.  and  whose  bases  are  identical  faces  in  the  planar  arrangements 
Qi^f)  and  Q(^c)  induced  in  A^  and  A,.,  respectively,  by  their  intersections  with  the 
vertical  rectangles.  Qi^f)  and  Q(^c)  each  contain  at  most  n  segments,  so  the  complex- 
ity of  a  single  face  in  either  of  the  two  arrangements  is  0(na{n)).  Since  C  is  a  prism, 
its  complexity  is  proportional  to  the  complexity  of  either  of  its  bases,  implying  our  first 
claim. 

There  are  two  types  of  horizontal  exposed  segments  on  the  boundaries  of  type  (2) 
cells:  portions  of  the  horizontal  edges  of  vertical  rectangles  and  portions  of  edges  of  the 
horizontal  plates.  We  will  restrict  our  attention  to  the  latter;  the  former  are  handled 
similarly.  Let  p  be  the  plane  containing  some  horizontal  plate  A.  Intersections  of  verti- 
cal rectangles  with  p  induce  a  planar  arrangement  Q  =  Q(A)  in  ;?  of  at  most  n  segments. 
Refine  the  type  (2)  cells  by  adding  a  wall  filling  each  face  of  Q  that  lies  outside  of  A 
and  is  incident  to  the  relative  boundary  3A  of  A.  Repeating  the  same  procedure  for 
every  horizontal  plate  (and  every  horizontal  edge  of  a  vertical  rectangle,  within  the  hor- 
izontal plane  containing  it),  we  obtain  a  decomposition  of  all  cells  of  type  (2)  into  sub- 
cells.  Note  that  all  horizontal  exposed  segments  have  been  eliminated,  so  each  subcell  is 
of  type  (1)  (and  thus  is  a  right  vertical  prism).  We  will  now  argue  that  the  bases  of  sub- 
cells  contained  in  type  (2)  cells  have  total  complexity  0(n^a{n)),  so  that  the  total  com- 
plexity of  these  subcells  is  bounded  by  0(n^a(n))  as  well;  this  would  complete  the 
proof  of  the  upper  bound  as  the  only  elements  of  type  (2)  cells  not  accounted  for  in  the 
combined  complexity  of  their  subcells  is  that  of  horizontal  exposed  segments  and  their 
incident  faces  and  vertices,  and  there  are  0(n^  +  q)  exposed  segments  altogether,  with 
each  exposed  segment  incident  to  exactly  two  faces  and  two  vertices. 

Consider  a  subcell  C  of  a  type  (2)  cell.  Let /and  c  be  the  floor  and  ceiling  bases 
of  C  contained  in  horizontal  planes  pf  and  pc,  respectively.  By  construction, 
dC  — /U  c  must  consist  of  portions  of  the  vertical  rectangles  extending  all  the  way  from 
/  to  c.  In  particular,  any  cross  section  of  C  by  a  horizontal  plane  /?  is  a  face  in  the 
planar  arrangement  Qip)  of  segments  induced  in  p  by  its  intersections  with  the  vertical 
rectangles.  Therefore  /  (respectively,  c)  is  the  union  of  one  or  more  such  faces  in  an 
analogous  arrangement  Q(pf)  in  pf  (respectively,  Q{pc)  in  pc),  which  is  identical  to 
Qip)  except  for  the  addition  of  the  edges  of  the  plate  contained  in  pf  (resp.  c). 
Observe  that,  by  construction,  the  closure  of  either  /or  c  (or  both)  must  contain  a  hor- 
izontal exposed  segment,  for  otherwise  C  would  coincide  with  an  original  cell  of  type 
(1).  Let  us  suppose  that  it  is  /  that  meets  some  horizontal  exposed  edge.  Thus  /coin- 
cides with  the  union  of  some  horizon  faces  of  the  boundary  3Ay  of  the  convex  plate  A^ 
in  the  arrangement  Qipf).  Moreover,  it  is  easily  checked  that  each  such  horizon  face  in 
QiPf)  lies  in  the  base  of  at  most  two  subcells  (in  the  floor  of  one  cell  and  in  the  ceiling 
of  another).  However,  the  results  of  [EGPPSS]  are  easily  seen  to  imply  that  the  horizon 
of  any  closed  convex  curve  in  a  planar  arrangement  of  n  segments  has  complexity 
0(na(n));  thus,  the  combined  complexity  of  the  bases  of  all  subcells  of  type  (2)  cells  is 
Oin^a{n)),  as  asserted. 

The  lower  bound  of  Cl(n  a(n)  +  q)  for  the  complexity  of  a  single  cell  is  easily 
obtained  by  a  slight  modification  of  the  construction  of  Theorem  2.5.  Thus  we  have 
shown 

Theorem  4.5:  In  an  arrangement  of  n  horizontal  convex  plates  and  vertical  rectangles, 
the  worst-case  complexity  of  a  single  cell  is  0(n^a(/i)  +  ^),  where  q  is  the  total  number 
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of  corners  of  the  horizontal  plates,   n 

5.   Algorithms 

Here  we  come  to  the  issue  of  actually  building  castles  in  the  air,  as  promised  in  the 
title  of  this  paper.  Namely,  we  are  interested  in  computing  one  or  more  cells  in  an 
arrangement  A  of  n  triangles  in  R^.  Note  that  it  is  not  difficult  to  compute  the  full 
arrangement  A  in  time  €K/i^)  by  modifying  the  algorithm  of  [EOS]  for  constructing  the 
arrangement  of  n  planes.  Thus  our  goal  is  to  attain  subcubic  performance  in  computing 
a  single  cell  (or  all  interesting  cells).  The  previous  results  of  [PS]  are  non-algorithmic, 
and  our  combinatorial  analysis  cannot  be  immediately  translated  into  an  efficient  algo- 
rithm either.  In  Section  5.1  we  describe  a  randomized  algorithm  for  computing  a  single 
cell  in  an  arbitrary  arrangement  of  n  triangles;  its  expected  time  complexity  is  asymptot- 
ically nearly  optimal.  Moreover,  the  same  algorithm  takes  only  0{n^'^^)  time,  for  any 
8>0  (with  the  constant  of  proportionality  depending  on  8),  when  applied  to  arrange- 
ments of  vertical  rectangles  and  horizontal  convex  plates  discussed  in  the  previous  sec- 
tion (assuming  that  each  horizontal  plate  has  constant  complexity),  or  to  arrangements 
of  n  triangles  lying  in  planes  with  a  constant  number  of  distinct  orientations.  We  then 
comment  on  the  difficulty  of  extending  our  algorithm  to  compute  many  cells,  and  out- 
line an  alternative,  less  efficient  but  still  subcubic  algorithm  for  obtaining  all  interesting 
cells,  which  is  based  on  the  recent  results  of  [EGHSSSW].  Finally,  we  present  a  custom- 
ized deterministic  algorithm  for  arrangements  of  triangles  with  few  distinct  plane  orien- 
tations (as  in  Section  4). 

All  of  our  algorithms  use  a  common  data  structure  to  represent  cells  in  three- 
dimensional  arrangements  of  triangles.  A  cell  is  represented  by  the  list  of  its  boundary 
components,  starting  with  the  explicitly  marked  outer  boundary  (unless  the  cell  is 
unbounded)  and  followed  by  zero  or  more  "islands",  all  linked  to  one  another  (the 
exact  nature  of  this  linkage  will  be  discussed  below).  Each  boundary  component  is 
stored  in  the  form  of  an  "adjacency  structure",  where  each  face  is  given  by  its  outer 
boundary  (unless  the  face  is  unbounded)  plus  zero  or  more  islands,  properly  linked  to 
one  another,  each  of  which  is  in  turn  represented  as  a  circular  list  of  incident  edges  and 
vertices.  Conceptually,  we  store  each  face  twice  (corresponding  to  the  two  sides  of  the 
triangle),  leading  to  similar  multiplicities  in  the  representation  of  edges  and  vertices,  as 
in  the  puffy  model  discussed  in  Section  2.1.  This  ensures  that  every  boundary  element 
is  incident  to  a  unique  cell  (duplicate  elements  are  still  linked  among  themselves;  in  fact 
an  actual  implementation  need  not  store  them  as  physically  distinct  entities).  Thus  each 
(copy  of  an)  edge  is  incident  to  exactly  two  faces,  a  triple-intersection  vertex  is  incident 
to  three  faces  and  three  edges,  etc.  The  precise  implementation  of  such  a  data  structure 
can  be  done  along  the  lines  suggested  in  [TDS]  or  in  [DL];  however  we  will  not  be  con- 
cerned here  with  such  implementation  details. 

5.1.   Calculating  a  single  cell 

Our  algorithm  proceeds  as  follows: 

Input:  A  set  G  oi  n  triangles  in  R^,  a  point  p,  and  an  "efficiency  parameter"  8>0. 

Output:  The  cell  of  A{G)  containing  p,  represented  by  the  incidence  structure  discussed 
above. 

(1)     Select  a  random  sample  R  of  r  triangles  (where  r  is  a  constant,  depending  on  8,  to 
be  determined  later). 
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(2)  Construct  the  full  arrangement  A(R)  (using  any  convenient,  even  brute-force, 
method). 

(3)  By  exhaustive  search,  identify  the  cell  Cr  of  A(R)  that  contains  p. 

(4)  Subdivide  Cr  into  open  tetrahedra  as  in  the  Slicing  Theorem  (again,  using  brute 
force) . 

(5)  Let  to  be  the  tetrahedron  containing  p.  For  each  triangle  A  ^  G  compute  the  convex 
polygon  A  n  to,  and  let  G,^  be  the  collection  of  such  (triangulated)  non-empty 
polygons.  Apply  the  algorithm  recursively  to  (G,^,p,8)  to  obtain  a  cell  Co-  If  Co  is 
a  bounded  cell  of  A(G,^),  stop  — C^  is  the  desired  cell  of  A(G)  containing  p.  Other- 
wise apply  the  trimming  step  (8)  below  to  C,^  =  Cg  (denoting  the  resulting  trimmed 
subcell  containing  p  by  Co)- 

Repeat  steps  (6)-(8)  for  each  of  the  remaining  tetrahedra  ti'to  in  the  decomposition  of 
Cr  and  proceed  to  step  (9). 

(6)  As  for  to,  compute  for  each  triangle  A  €  G  the  convex  polygon  A  D  r,  and  let  G,  be 
the  collection  of  all  such  non-empty  (triangulated)  polygons. 

(7)  Apply  the  algorithm  recursively  to  (G,,poo,8)  where  p«  is  'the  point  at  infinity',  to 
obtain  the  unbounded  cell  C,  of  A(G,). 

(8)  Trim  C,  along  the  faces  of  t  and  discard  the  outer  portion  of  the  cell;  this  may  split 
C,  into  several  subcells  all  lying  within  t  and  incident  to  the  boundary  dt  of  t. 

A 

(9)  Reconstruct  the  cell  C  of  A(G)  containing  p  by  starting  with  Co  and  'growing'  it  as 
follows:  Let  a  window  of  a  (trimmed)  cell  C  of  A  (G,)  be  a  face  of  C  on  dt  intro- 
duced by  the  trimming  step  (8).  Clearly  each  window  is  shared  by  two  cells  lying 
in  adjacent  tetrahedra^  A  window  is  transparent  if  it  is  not  contained  in  a  triangle 
of  G.  Starting  with  Co,  locate  all  cells  sharing  transparent  windows  with  it  and 
'glue'  them  to  Co  along  the  windows.  Repeat  until  the  resulting  cell  C  has  no 
unpaired  transparent  windows.   This  is  the  desired  cell. 

Why  does  the  algorithm  work?  Observe  that  the  cell  Cq  of  A(G)  containing  p  is 
completely  contained  in  the  cell  Cr  of  A{R)  containing  p,  so  that  the  algorithm  correctly 
restricts  its  computation  to  Cr.  Notice  that  if  p  does  not  lie  in  the  unbounded  cell  C,^  of 
A(Gfg)  then  Cq  must  coincide  with  Co,  in  which  case  no  further  computation  is  neces- 
sary and  the  algorithm  correctly  returns  Co-  Assume  p  €Ct^  =  Co  which  is  unbounded. 
If  a  point  X  lies  in  Cq  ~  Co,  there  is  a  path  ct  in  Cq  connecting  it  to  p.  A  traversal  of  ct 
from  p  to  X  starts  off  in  Co  and  proceeds  to  visit  trimmed  cells  of  A(G,)  (for  various  r), 
moving  from  cell  to  cell  through  transparent  windows.  Moreover,  each  cell  C  crossed 
by  a  must  be  a  portion  (within  the  corresponding  tetrahedron  r)  of  the  unbounded  cell 
C,  of  A(G(),  because  C  can  be  reached  from  outside  r.  Hence  x  will  be  in  the  cell  built 
by  the  algorithm.  So  the  constructed  cell  is  no  smaller  than  Co  and.  trivially,  it  cannot 
be  larger  than  Cc  either.   Thus  the  cell  the  algorithm  computes  is  precisely  Cq- 

We  next  analyze  the  expected  running  time  of  the  algorithm.  Since  r  is  a  constant. 
all  steps  of  the  algorithm  besides  recursive  invocations  and  the  trimming  and  reconstruc- 
tion steps  can  be  performed  in  overall  linear  time.  Denoting  by  \C  \  the  complexity  of  a 
cell  C,  observe  that  trimming  the  unbounded  cell  C,  of  AiG,),  for  some  tetrahedron  t, 
can  be  performed  in  time  proportional  to  \C,\,  as  it  is  sufficient  to  check  each  feature  of 
C,  against  each  face  of  t.    The  time  required  for  the  reconstruction  step  is  proportional 
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to  the  total  complexity  of  all  transparent  windows,  because  this  step  can  be  accom- 
plished by  simultaneously  tracing  the  common  face  of  each  pair  of  adjacent  tetrahedra 
using  the  incidence  structure  of  the  trimmed  cells.  Determining  whether  Co  is  the 
unbounded  cell  of  A(G,^)  can  be  accomplished  in  constant  time,  since  the  representation 
we  use  explicitly  stores  the  outer  component  (if  any)  of  a  cell  boundary.  Hence  the 
total  time  complexity  of  all  steps  of  the  algorithm,  excluding  recursive  calls,  is  bounded 

by 

t 

for  some  positive  constants  A,  B  depending  on  r,  with  the  summation  taken  over  all 

tetrahedra  t  contained  in  Cr.    Observe  that  a  small  perturbation  of  UG,  (that  moves 

r 

tetrahedra  away  from  each  other)  produces  an  arrangement  of  2  l^tl  triangles  in  which 

t 

the     boundary     of     the     unbounded     cell     is     precisely      UdC,.       In     particular, 

2  |C,  I  ^^(2  l<^»|)-    It  immediately  follows  that  the  time  complexity  of  the  algorithm 
r  t 

exclusive  of  recursive  calls  is  bounded  above  by  0(/i  +  iC^  \G,\)).   Now  ^  \G,\  is  easily 

t  t 

seen  to  be  bounded  by  Dn,  for  some  constant  D  depending  on  r.  Therefore,  the  algo- 
rithm requires  0(C(Dn))  overhead  and  recurs  on  K  problems  (G,,,pco,8),  where  K  is  the 
total  number  of  tetrahedra  contained  in  Cr  (by  the  Slicing  Theorem,  K^Jl,ir)  where  7 
is  a  constant  independent  of  r),  and  /,-  varies  over  all  such  tetrahedra  (including  to).  By 
the  €-net  theory  of  [HW]  (see  also  [CI]),  with  high  probability,  each  tetrahedron  meets 

at  most  logr  triangles,  for  some  constant  a>0  independent  of  r.    Moreover,  with 

no  extra  overhead  we  can  verify  that  this  does  indeed  hold  for  our  sample  /?;  if  not  we 
simply  discard  R  and  try  another  sample,  until  we  hit  one  which  has  the  €-net  property. 
It  easily  follows  that  the  expected  number  of  such  iterations  is  a  constant  (depending  on 
r).     Hence,    in   expected   linear   time,   the   algorithm   produces   a   sample   for   which 

\G,i\  ^ logr,  for  all  i.   In  particular,  denoting  by  T(n)  the  expected  running  time  of 

the  algorithm,  we  obtain  the  following  recurrence  for  T(n): 

T{n):< 


Kn^\ogr)  +  Oa(Dn)),  if  n>r. 


otherwise. 
Define  7  =  lim  sup      °  ^\  '  .    Recall  that  i(n)  has  been  shown  to  be  Q,(n^a(n))  and 

h^co  lOgrt 

0(n'^'^a(n)'^'^\og'^'^n),  so  2  ^  7  <  7/3.  Hence,  for  any  8  >  0,  there  is  a  choice  of  a  con- 
stant b  such  that  the  term  0(t,(Dn))  is  bounded  above  by  bn^'^^'^  for  all  n  (with  b 
depending  on  8  and  r).   Thus  the  recurrence  reduces  to 


r(n)< 


an 

r 


KT(^\ogr)+bn^^^'^,  for  «  >  r, 

otherwise, 


which  has  a  solution  T{n)<Hn°-  for  any  a>max(7  +  8/2,  logAT/log— 7^ — )  (with  H 

a  logr 
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depending  on  r,  8  and  a).  Notice  that  the  second  quantity  can  be  made  not  to  exceed 
7  +  8/2  by  choosing  a  sufficiently  large  r,  as  K^Jt,{r).  Choosing  such  an  r  for  the 
given  input  value  of  8,  we  obtain  T(n)  =  0(n^'^^),  with  the  constant  depending  on  8. 
In  particular,  we  have  shown: 

Theorem  5.1:  Given  a  fixed  8>0,  a  single  cell  in  an  arrangement  of  n  triangles  in  R^ 
can  be  computed  in  expected  time  Cj/i"'"^*,  with  the  constant  Cg  depending  on  8,  where 

"y  =  limsup  '°^^W  ^[2.7/3].   □ 
t-=o        log/: 

Corollary  5.2:  A  single  cell  in  an  arrangement  of  n  triangles  can  be  computed  in 
expected  time  0(n^''^'^^),  for  any  8>0,  with  the  constant  of  proportionality  depending 
on  8. 

Notice  that,  should  the  conjecture  posed  in  the  introduction  be  demonstrated  to 
hold,  the  above  theorem  would  immediately  provide  an  0(_n^'^^)  expected  time  algo- 
rithm for  computing  any  single  cell  of  the  arrangement. 

We  also  note  that  the  polygons  passed  to  recursive  invocations  of  the  algorithm  will 
not  in  general  be  triangles.  However,  it  is  easy  to  see  that,  independently  of  the  level 
of  recursion  on  which  such  a  polygon  is  created,  it  can  be  represented  as  the  intersection 
of  an  original  triangle  with  a  single  tetrahedron.  Hence  such  a  polygon  can  have  no 
more  than  7  sides,  and  the  presence  of  non-triangles  can  be  compensated  for  by  using 
i{5n)  instead  of  C(i)  in  our  argument  (as  any  convex  polygon  with  up  to  7  sides  can  be 
cut  into  at  most  5  triangles),  and  by  replacing  the  constant  a  by  5a. 

Remark:  Several  recent  techniques  (due  to  Clarkson  [CI],  Chazelle  and  Friedman  [CF], 
and  Matousek  [Ma])  provide,  in  certain  special  cases,  tools  either  for  a  deterministic 
construction  of  a  triangulation  of  space  with  properties  similar  to  those  discussed  above, 

or  for  obtaining  such  a  decomposition  in  which  each  cell  is  cut  by  only  0(— )  objects, 

rather  than  0(— log  r).    We  do  not  know  whether  these  techniques  can  be  adapted  to 

the  problem  at  hand,  so  as  to  make  the  above  algorithm  deterministic  or  further 
improve  its  time  complexity. 

Finally,  observe  that  the  time  complexity  analysis  of  the  above  algorithm  relies  on 
the  following  facts: 

(1)  The  complexity  of  any  cell  in  the  arrangement  is  bounded  by  a  function  (,(n)  such 
that  a  cell  in  the  arrangement  formed  by  any  subset  of  r  triangles  (a)  has  complex- 
ity 4(r)  and  (b)  can  be  triangulated  into  0(t,ir))  tetrahedra. 

(2)  With  high  probability,  the  number  of  triangles  cutting  any   tetrahedron   that   is 

missed  by  the  r  sample  triangles  is  at  most  logr. 

r 

(3)  The  intersection  of  a  tetrahedron  with  a  triangle  is  a  convex  polygon  of  constant 
complexity. 

Thus,  in  many  situations  where  the  above  conditions  hold  (perhaps  with  a  different 
function  vj;  taking  the  place  of  Q,  the  algorithm  will  correctly  compute  the  desired  cells 
in  expected  time  0{i\i(n)n^),  for  any  8>0.    In  particular,  we  have: 

Theorem  5.3:  Given  a  fixed  8  >  0,  a  single  cell  in  an  arrangement  of  n  triangles  lying 
in  planes  with  a  constant  number  /of  distinct  orientations  can  be  computed  in  expected 
time  Cg./n        ,  with  the  constant  Cj,/ depending  on  8  and  /. 
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Theorem  5.4:  Given  a  fixed  8  >  0,  a  single  cell  in  an  arrangement  of  n  horizontal  trian- 
gles and  vertical  rectangles  can  be  computed  in  expected  time  CgAi^"*"^,  with  the  constant 
Cg  depending  on  8. 

Observe  that  Theorem  5.3  is  superceded  by  Theorem  5.6  of  Section  5.3  which,  for 
a  constant  /,  provides  a  faster  and  more  general  deterministic  algorithm  for  computing 
portions  of  such  arrangements. 

5.2.   Calculating  many  cells 

An  attempt  to  extend  the  algorithm  of  Section  5.1  to  calculate  an  arbitrary  collec- 
tion of  cells  of  an  arrangement  of  triangles  (e.g.  all  interesting  cells)  faces  the  following 
technical  problem:  our  algorithm  takes  advantage  of  the  fact  that  the  desired  single  cell 
is  contained  in  a  single  cell  of  the  arrangement  A(i?)  of  the  random  sample  of  triangles. 
This  implies  that  the  number  of  tetrahedra  that  require  further  processing  is  only 
Oir'''^air)^'^\og^'^r)  (in  fact,  0(t,(r))),  which  leads  to  a  recurrence  with  a  favorable 
time  complexity.  In  contrast,  when  calculating  many  cells,  there  is  no  a  priori  sharp 
bound  on  the  overall  complexity  of  the  cells  of  A  (/?)  which  contain  the  desired  cells  of 
A(G)  — in  the  worst  case  all  cells  of  A(/?)  could  be  involved— so  that  no  comparable 
recurrence  relation  can  be  obtained.  Hence,  although  our  algorithm  can  be  easily 
adapted  to  calculate  a  collection  of  cells,  we  do  not  have  sharp  worst-case  bounds  on  its 
expected  running  time  in  this  case. 

Instead,  we  can  use  the  following  alternative  approach.  Edelsbrunner  et  al. 
[EGHSSSW]  recently  described  a  randomized  procedure  which,  given  a  planar  arrange- 
ment Q  of  n  segments,  preprocesses  it  in  0(n^'^\og^'^n)  expected  time  and 
0(/i'*'^a(n)log/i)  space  so  that,  given  any  query  point  x,  one  can  calculate  the  face  of  Q 
containing  x  in  time  0(n^'^log^^'^n  +  k),  where  k  is  the  complexity  of  that  face.  Using 
this  procedure,  a  subset  of  the  cells  of  A(G)  can  be  computed  as  follows. 

Consider  first  the  task  of  calculating  all  interesting  cells.  We  begin  by  obtaining  a 
point  on  each  exposed  segment  (in  overall  O(n^)  time).  Consider  each  of  these  points  y 
in  turn,  and  let  A  be  a  triangle  containing  y.  The  algorithm  proceeds  to  trace  dC,  where 
C  is  the  cell  containing  y,  as  follows.  It  first  obtains  the  face  /  containing  y  in  the 
arrangement  induced  in  A  by  intersections  with  the  remaining  triangles,  using  the  pro- 
cedure of  [EGHSSSW],  and  then  moves  to  adjacent  faces  of  dC,  which  share  an  edge 
with  /and  lie  in  other  triangles,  obtains  each  of  them,  and  repeats  this  procedure  until 
the  entire  boundary  component  is  traced  this  way.  One  then  needs  to  'hop'  from  one 
boundary  component  of  C  to  another.  This  is  done  by  'linking'  triangles  in  the  follow- 
ing fashion:  In  the  preprocessing  stage,  for  each  of  the  3n  triangle  corners,  one  com- 
putes (by  brute  force)  the  triangle  that  lies  directly  below  it  (if  any)  and  the  point  of 
intersection  of  this  triangle  with  a  vertical  line  through  the  corner.  Given  such  a  point  z 
lying  directly  below  a  corner  c,  one  proceeds  (again,  by  brute  force)  to  the  boundary  of 
the  face  containing  z  (in  the  induced  planar  arrangement)  and  then  follows  the  boundary 
until  the  first  vertex  c'  of  that  face  is  encountered;  then  c'  is  linked  to  c  and  vice  versa. 
Corners  that  have  no  triangle  below  them  are  linked  to  the  'vertex  at  infinity'.  The 
desired  "boundary-hopping"  is  then  accomplished  by  examining  all  vertices  in  the  com- 
puted boundary  component  and  following  the  links  from  any  of  them  to  vertices  that 
have  not  yet  been  encountered,  which  then  serve  as  starting  points  y  in  new  boundary 
components.  Consider  an  interior  boundary  component  AT  of  a  bounded  cell.  Since  the 
lowest  vertex  of  ^  is  a  triangle  corner,  it  is  easy  to  verify  that  some  link  connects  K  to 
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either  the  outer  boundary  component  or  to  an  interior  component  whose  lowest  point  is 
below  that  of  K.  This  immediately  implies  that  all  interior  components  are  (possibly 
indirectly)  linked  to  the  outer  component.  This,  in  particular,  allows  one  to  pass  from 
one  boundary  component  of  a  cell  to  another  until  all  have  been  traced.  The  argument 
in  the  case  of  the  unbounded  cell  is  similar,  with  the  'vertex  at  infinity'  playing  the  role 
of  the  outer  component  of  the  boundary.  It  is  easy  to  verify  that  linear  time  is  sufficient 
to  compute  each  of  the  links  by  brute  force,  thus  only  quadratic  additional  time  suffices 
to  guarantee  that  no  component  of  the  desired  boundary  will  be  left  unvisited. 

Consider  next  the  general  case  where  we  want  to  calculate  an  arbitrary  collection  of 
cells,  each  specified  by  a  point  x  in  its  interior.  For  each  such  cell  C,  we  first  identify 
some  point  y  on  its  boundary  dC.  This  can  be  done  in  linear  time  per  cell  by  choosing  y 
to  be  the  first  point  of  intersection  of  the  downward-directed  vertical  ray  emanating 
from  X  with  any  triangle;  if  there  is  no  such  point,  let  y  be  the  'vertex  at  infinity'-C  in 
this  case  is  the  unbounded  cell  and  the  algorithm  will  then  start  at  the  vertex  at  infinity 
and  follow  all  precomputed  links  from  it  to  boundary  components  of  C  (as  in  the  discus- 
sion just  presented). 

If  the  number  m  of  points  marking  cells  is  very  large,  this  linear  overhead  per 
marker  x  may  be  too  expensive.  In  this  case  we  can  locate  the  desired  points  y  as  fol- 
lows. Extend  each  triangle  into  a  full  plane,  and  apply  the  (randomized)  algorithm  of 
[EGS3]  which  calculates  the  plane  lying  immediately  below  each  of  these  m  points  in 
expected  time  0{m^'^~^n^''^  \o^n)  for  any  8>0.  We  thus  obtain  m  points  y,  each 
lying  below  a  corresponding  given  point  x  on  one  of  the  n  planes  containing  the  trian- 
gles. Using  the  precomputed  data  structures  of  [EGHSSSW]  for  fast  face  queries,  we 
compute,  for  each  y,  a  point  on  the  boundary  of  the  face  /containing  y  in  the  arrange- 
ment induced  in  the  plane  of  the  corresponding  triangle,  thereby  obtaining  a  point  on 
the  boundary  of  the  cell  containing  y  (and,  therefore,  x).  This  requires  0{n^''^\og^^'^n) 
expected  time  per  query. 

Having  obtained  in  this  manner  an  "anchor  vertex"  on  the  boundary  of  each  of  the 
desired  cells,  we  can  then  continue  as  in  the  case  of  all  interesting  cells  described  above. 
We  omit  further  details  concerning  this  extension.  Thus  we  obtain 

Theorem  5.5:  Any  collection  of  cells  of  A{G),  defined  by  specifying  a  set  of  m  points 
marking  the  desired  cells,  can  be  computed  in  expected  time 

0(m3''*-8n3'^^3»log2«  -H  (A/  +  mW^og^^'^n  ^  n^'^og^'^n) 

and  0(m3"*-*n^"''^38  ^- m  +  n"^ '"^ a{n)\ogn)  working  storage,  for  any  8>0  (with  the  con- 
stant of  proportionality  depending  on  8),  where  M  is  the  total  complexity  of  the  desired 
cells.  If  the  number  m  of  marking  points  is  small,  a  more  straight  forward  approach 
will  compute  the  marked  cells  in  expected  time 

0{mn  +A/n^'3log"'2„  ^  n^'^\og^''^n) 

using  O {M  +  n^ '^ a{n)\ogn)  space.  A  somewhat  simpler  algorithm  computes  all  interest- 
ing cells  of /\(G)  in  expected  time  0(,n^'^OL{n)'^'^\og^^'^n)  and  space  0(n^'^a(n)logn). 

Remark:  Notice  that  the  amount  of  work  required  for  computing  an  arbitrary  subset  of 
cells  is  largely  dominated  by  the  preprocessing  time,  which  in  turn  requires  nearly  as 
much  time  as  that  for  computing  all  interesting  cells.  Therefore,  it  may  be  assumed  that 
all  interesting  cells  are  to  be  calculated  anyway  and  one  is  facing  the  problem  of  com- 
puting some  additional  dull  (i.e.,  convex)  cells  of  an  arrangement  of  triangles,  where 
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each  cell  is  specified  by  a  point  in  it. 

5.3.   A  detenmnistic  algorithm  for  arrangements  with  few  orientations 

In  this  section  we  sketch  an  efficient  0((/n  +  M  +  «^)/logn)-time  deterministic 
algorithm  for  computing  a  subset  of  m  cells  in  an  arrangement  of  n  convex  plates  lying 
in  planes  with  only  /distinct  orientations,  where  M  is  the  total  complexity  of  the  cells 
being  computed.  For  simplicity  of  presentation  we  will  assume  that  the  plates  are  in 
fact  triangles  and  that  no  more  than  three  triangles  meet  at  a  common  point,  but  a  more 
general  situation  can  be  handled  just  as  easily.  Each  desired  cell  C  is  identified  by  a 
point  xiC.  Given  x,  the  algorithm  locates  a  vertex  v^  on  dC,  using  a  technique  similar 
to  that  in  Section  5.2,  but  admitting  a  much  simpler  implementation.  Specifically,  the 
plane  defined  by  a  triangle  that  lies  immediately  below  x  can  be  determined  in 
Oiflogn)  time  by  performing  /  binary  searches,  one  in  each  collection  of  parallel 
planes.  Once  the  point  y  lying  in  such  a  plane  immediately  below  x  is  located,  a  point  on 
dC  is  easily  obtained  by  one  O (/log n) -time  ray-shooting  query  (described  below)  in 
that  plane,  and  a  vertex  on  dC  can  be  obtained  after  two  more  queries  of  this  kind). 
Then  it  traces  out  all  edges  and  vertices  of  dC  by  starting  at  Vg  and  repeatedly  "shoot- 
ing" along  an  edge  of  dC  incident  to  the  "current"  vertex  v  and  thus  discovering  a  new 
vertex  w  of  dC.  This  gives  us  two  new  edges  incident  to  w  and  we  shoot  along  them  to 
discover  further  vertices. 

Two  specific  problems  arise  in  this  approach.  First,  how  can  each  "shot"  be  per- 
formed in  0(f\ogn)  time?  Shooting  along  exposed  segments  is  easily  implemented  by 
precomputing  the  answers  to  all  possible  queries,  while  each  ray-shooting  query  along  a 
line  of  intersection  of  two  plates  reduces  to  a  two-dimensional  ray-shooting  in  some 
arrangement  induced  in  either  plate  by  its  intersections  with  the  remaining  plates.  By 
assumption,  each  such  arrangement  A,-  in  plate  A,-  is  formed  by  at  most  /—I  overlaid 
subarrangements,  each  of  which  consists  of  a  collection  of  parallel  segments.  Moreover, 
there  are  only  /—I  possible  direction  for  the  shooting  ray.  Thus,  for  each  subarrange- 
ment  and  each  shooting  direction  we  can  prepare  a  data  structure  that  supports  0(log/j) 
shooting  queries  using,  for  example,  the  technique  of  [ST].  It  is  easily  checked  that  this 
approach  guarantees  0(flogn)  query  time  while  the  preprocessing  can  be  accomplished 
for  all  n  plates  in  0(fnhogn)  time  and  0(Jh^)  storage.  (This  also  takes  care  of  the  ini- 
tial step  of  moving  from  each  given  point  x  to  a  vertex  on  the  boundary  of  the  cell  con- 
taining it.) 

The  second  problem  that  arises  is  handling  clusters  of  edges  lying  in  dC  but  not 
connected  to  each  other.  Such  situations  can  occur  if  there  are  'islands'  in  faces  of  dC, 
or  if  dC  is  not  connected.  We  have  dealt  with  the  latter  case  in  Section  5.2  by  linking 
each  triangle  corner  to  (a  vertex  lying  in)  the  triangle  directly  below  it  to  facilitate  "hop- 
ping" from  one  component  of  dC  to  another.  The  former  situation  is  approached  simi- 
larly except  that  the  linking  is  performed  in  each  planar  arrangement  A,-  induced  in  the 
plate  A,-  by  intersections  with  the  remaining  plates.  For  each  segment  endpoint  c  in  A,-, 
we  locate  the  segment  s  lying  immediately  below  c  (in  A,).  Then  c  is  linked  to  the  ver- 
tex of  Ai  lying  on  s  nearest  to  the  point  of  s  immediately  below  c.  If  the  "linking"  steps 
are  carried  out  using  ray-shooting,  the  total  overhead  can  be  shown  not  to  exceed 
OifiiHogn). 

Theorem  5.6:  Given  a  set  of  m  marking  points,  the  cells  marked  by  these  points  in  an 
arrangement  of  n  triangles  lying  in  planes  with  /distinct  orientations  can  be  computed  in 
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deterministic  time  0((m  +  M  +  n^)f\ogn)  and  0(M  +fn^)  storage,  where  M  is  the  total 
complexity  of  the  marked  cells.  In  particular,  if  no  two  markers  lie  in  the  same  cell,  the 
time  complexity  of  the  algorithm  reduces  to  0((M  +  n^)/logn).   □ 

6.   Discussion  and  open  problems 

The  paper  leaves  one  major  open  problem,  namely  the  settling  of  our  conjecture 
that  5(n)  =  €Kn^a(n)).  In  addition,  we  note  the  following  extensions  and  applications 
of  our  results: 

(1)  As  already  hinted,  the  analysis  of  Section  2  easily  extends  to  yield  a  similar  upper 
bound  for  the  complexity  of  all  non-convex  cells  in  an  arrangement  of  n  arbitrary  con- 
vex (flat)  plates  in  3-space  (with  an  additive  correction  term  that  accounts  for  the  com- 
plexity of  the  plate  boundaries).   We  omit  the  straightforward  details  of  this  extension. 

(2)  We  believe  that  our  results  can  be  generalized  to  higher  dimensions;  specifically  we 
conjecture  that  the  maximum  number  of  facets  (highest-dimensional  faces)  of  all  non- 
convex  cells  in  an  arrangement  of  n  ^-simplices  in  (^f  +  l)-dimensional  space  is  close  to 
0(n'^'^^'^).  (Note  that  the  lower  bound  construction  can  be  generalized  to  any  number 
of  dimensions,  yielding  a  lower  bound  of  ^(n''''"^'^)  for  this  complexity.)  The  technique 
used  in  [EGS3]  to  obtain  such  an  extension  for  the  case  of  hyperplanes  may  be  useful 
for  proving  this  conjecture. 

(3)  The  technique  we  used  to  prove  the  Combination  Lemma  can  be  adapted  to  yield 
simpler  proofs  for  other  combination  lemmas,  such  as  those  given  in  [EGS2],  [EGS3]. 
Appendix  B  exemplifies  this  claim  by  providing  a  simple  proof  for  a  variant  of  the  two- 
dimensional  combination  lemma  of  [EGS2]  for  faces  in  arrangements  of  segments. 

(4)  Expanding  upon  our  motion-planning  application,  we  note  that  our  analysis  shows 
that  the  accessible  portion  C  of  the  configuration  space  FP  of  a  general,  not  necessarily 
convex,  polyhedron  B  with  K  faces  translating  amidst  polyhedral  obstacles  having  n 
faces  altogether,  has  0({Kn)^''^a(Kn)^''^\og*'^{Kn))  combinatorial  complexity,  and  that  it 
can  be  calculated  in  expected  time  O(iKn)^''^'^^)  (more  precisely,  in  time 
0{i(Kn)iKn)^)),  for  any  8>0.  However,  to  actually  plan  a  motion  between  two  given 
placements  zi,  Z2  of  fl,  it  may  not  even  be  necessary  to  calculate  the  entire  component  C 
of  FP  containing  these  placements.  In  this  regard  we  claim  that  if  B  can  translate  from 
zi  to  Z2,  then  it  can  always  do  it  along  a  polygonal  path  having  only  0{{Kn)^)  turns, 
and  that  in  the  worst  case  (at  least  for  a  constant  K)  that  many  turns  are  necessary.  The 
claim  immediately  follows  from  the  fact  that  C  can  be  cut  into  0{(,Kn)^)  convex  polyhe- 
dra,  as  shown  in  Section  3  (Corollary  3.2  to  the  Slicing  Theorem).  An  example  where 
£)iKn^)  turns  are  necessary  (for  some  non-convex  B)  is  not  difficult  to  construct.  How- 
ever, we  do  not  know  how  to  calculate  such  a  path  efficiently,  without  obtaining  first 
the  entire  component  C.  Further  extensions  and  applications  of  the  combinatorial 
results  obtained  above  to  motion  planning  problems  are  discussed  in  a  forthcoming  com- 
panion paper  [AS]. 

(5)  An  interesting  by-product  of  the  Slicing  Theorem  and  the  preceding  remark  is  that 
the  techniques  of  [SS]  for  coordinated  motion  planning  for  two  independent  robots  can 
be  applied  to  obtain  an  0{n^'^'^^^)  (actually,  0((C(n))^n^))  algorithm  (for  any  5>0) 
for  coordinating  the  motions  of  two  independent  translating  polyhedra,  each  having 
0(1)  complexity,  amidst  stationary  polyhedral  obstacles  with  a  total  of  n  faces.  See  [SS] 
for  more  details. 
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Appendix  A.  The  topology  of  arrangements  of  triangles 

In  this  appendix  we  derive  several  basic  properties  of  the  topological  structure  of 
an  arrangement  of  n  triangles  in  R^.  We  begin  with  the  analysis  of  the  number  of 
cutting-but-not-splitting  faces,  required  in  the  proof  of  the  Combination  Lemma  and 
proceed  to  justify  our  practice  of  bounding  the  complexity  of  a  cell  by  the  number  of  its 
faces.  Our  arguments  make  use  of  the  "puffy"  model  of  the  arrangement,  mentioned  in 
Section  2,  and  further  discussed  below. 

A.l  The  genus  of  cell  boundaries 

In  this  section  we  will  analyze  the  topological  structure  of  the  union  of  a  family  of 
triangles  in  R-'  in  order  to  bound  the  number  of  'cutting-but-not-splitting'  faces  in  the 
proof  of  the  Combination  Lemma. 

For  the  following  analysis,  it  will  be  useful  to  regard  the  boundary  of  a  cell  in  an 
arrangement  of  triangles  in  R''  as  a  compact  orientable  two-manifold.  To  this  end,  we 
replace  each  triangle  A  in  the  arrangement  by  a  "puffy  triangle"  A*  which  is  a  thin 
nearly  flat  body  bounded  by  two  surfaces,  one  on  each  side  of  A,  slightly  deformed 
away  from  one  another.  A*  has  the  same  edges  and  vertices  as  A,  and  is  sufficiently 
thin  so  that  the  combinatorial  pattern  of  intersections  of  the  puffy  triangles  is  identical 
to  that  of  the  original  arrangement  (such  choice  of  thicknesses  is  always  possible  under 
the  general  position  assumption).  Cells  in  the  new  arrangement  are  connected  com- 
ponents of  the  common  exterior  (i.e.,  complement  of  the  union)  of  all  puffy  bodies. 
However,  it  is  easy  to  check  that  now  cell  boundaries  are  compact  orientable  two- 
manifolds  (without  boundary).    Note  that  in  this  puffy  model,  every  original  face  occurs 
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in  the  new  arrangement  exactly  twice,  each  (pairwise  intersection)  edge— four  times, 
each  (triple-intersection)  vertex -eight  times,  and  each  vertex  that  is  the  intersection  of 
an  exposed  edge  with  another  triangle -twice.  Triangle  corners  and  exposed  segments 
are  not  duplicated. 

Let  us  recall  some  facts  from  elementary  algebraic  topology  (see  [GP]  for  refer- 
ence). For  a  cell  complex  X,  let  x(^)  be  its  Euler  characteristic.  It  has  the  property 
that,  for  any  two  cell  complexes  X  and  Y 

X(XUY)  =  xiX)  +  x(Y)-x(XnY)  (*) 

(where  X  D  y  is  a  sub-complex  of  X  and  Y).  We  will  require  the  following  standard 
facts:  x(-5^)=0,  x(S^  ^  0=0,  x('5^)  =  2,  x(0^)  =  l.  where  /=[0,1]  is  the  unit  interval, 
5^-a  circle,  5^-a  sphere,  D^-a  closed  disk,  and  S^  x  /-a  cylinder  (more  precisely,  the 
surface  of  a  bounded  cylinder  without  the  "lids").  In  particular,  notice  that  Euler 
characteristic  is  additive  when  X  and  Y  are  disjoint  or  intersect  along  a  set  of  disjoint  cir- 
cles and  and/or  cylinders.  Also  recall  that  the  Euler  characteristic  of  a  compact  con- 
nected orientable  two-manifold  (without  boundary)  is  equal  to  2(1  —  ^),  where  g  is  its 
genus. 

We  now  return  to  arrangements  of  triangles.  Consider  a  triangle  A  cutting  an 
arrangement  A  =A(G)  of  n  other  triangles.  The  intersection  of  A  with  A  is  a  planar 
arrangement  Q  =  Q(A)  of  (at  most)  n  segments  all  contained  within  A.  Keeping  in 
mind  the  "puffy"  model  of  the  arrangement,  observe  that  each  face  f  of  Q  corresponds 
to  two  faces  in  the  arrangement  A'  =  A(G  U{A}),  and  that  the  introduction  of  this  pair 
of  faces  changes  some  cell  C  of  A  (either  by  splitting  it  into  two  subcells,  or  by  just 
modifying  its  boundary  dC).  Our  goal  is  to  measure  the  effect  of  adding /on  the  genus 
of  dC.  For  the  remainder  of  this  discussion,  we  will  assume  that  A  does  intersect  some 
triangle  of  G,  for  otherwise  it  simply  adds  an  extra  component  (of  genus  0)  to  dC 
without  changing  the  rest  of  dC.  Let  us  denote  by  g  (respectively,  g')  the  genus 
summed  over  all  connected  components  of  dC  (respectively,  dC ,  where  C  is  the  one  or 
two  cells  produced  by  cutting  C  with  f),  and  by  c  (respectively,  c')  the  number  of  such 
connected  components.  As  we  will  only  be  interested  in  the  differences  8c  =  c'  —  c  and 
8^  =  g'  —  g,  we  may  assume  that  the  quantities  g,  c,  g' ,  and  c'  refer  only  to  the  com- 
ponents of  dC  and  dC  met  by  /.  Also  notice  that,  modulo  the  above  assumption,  the 
Euler  characteristic  of  dC  can  be  computed  as  xi^C)  =  2{c  —  g)  (simply  by  summing  x 
over  all  components  of  dC).  component  of  dC  separately  and  observing  that  the  Euler 
characteristic  is  additive  here  as  the  components  are  disjoint  by  definition).  A  similar 
identity  holds  for  x(5C")-   Two  cases  are  possible:  either 

(i)     /touches  the  relative  boundary  dA.  of  A,  in  which  case  the  introduction  of /cannot 
split  C  into  two  subcells,  or 

(ii)    /lies  within  the  relative  interior  of  A,  in  which  case  C  may  or  may  not  be  split  by  / 
into  two  disjoint  cells. 

Recall  that  the  boundary  of  /  consists  of  a  single  outer  component  and  zero  or  more 
islands,  all  enclosed  by  the  outer  boundary  and  enclosing  mutually  disjoint  regions  (that 
lie  outside  of/). 

We  will  consider  case  (i)  first.  Let  the  outer  boundary  of /be  denoted  douj.  Sup- 
pose that  /has  i^O  islands  and  that  doutf~dA  consists  of  b^O  connected  components. 
Since  we  have  assumed  that  A  intersects  some  triangles  of  G,  at  least  one  of  i,  b  is  not 
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zero.  In  the  puffy  model,  /can  be  identified  with  a  sphere  with  2i+b  open  discs 
deleted  (two  discs  are  deleted  for  each  island -one  on  either  side  of /—and  one  disk  for 
every  connected  component  of  doutf~^^)-   Using  (*),  we  obtain 

XC0  =  X(5^)  -  (2z  ^b)x{D'^)  +  (2/  +fe)x(-S^)  =  2  -  (2:  +b). 

Now  observe  that  the  boundary  dC  of  the  new  cell  C  can  be  obtained  by  removing  b 
disks  (one  for  each  contact  of  douJ  with  dC)  and  i  cylinders  (one  for  each  island  of  f) 
from  dC  and  "gluing"  the  remaining  portion  to  /along  the  li+b  "seam"  circles.  The 
Euler  characteristic  of  dC  with  b  disks  and  /  cylinders  removed  is  easily  seen  to  be 

XOC)-zx(5ix/)-fox(D^)  +  (2/+&)x(5^)  =  xOC)-Z>. 

Applying  (*)  once  again,  we  obtain  the  Euler  characteristic  of  dC: 

X(dC')  =  (x(dC)-b)  +  x(f)-(2i+b)xiS^)  =  X(dC)  +  2(\-i-b). 

Recalling  the  relation  between  the  genus  and  the  Euler  characteristic  of  a  compact 
orientable  two-manifold,  we  conclude  2(c'  —  g')  =  2ic  —  g)  +  2(1  —  i  —  b),  or 
8^  =  he  +  i  +  b  —  \.  Since  c'  (the  number  of  components  of  dC  met  by  /)  is  1  and 
1  <  c  <  I  +  fe,  it  follows  that  1  -  i  -  Z>  <  5c  <  0  and  thus  0  <  hg  <  i  +  b  -  I. 

Let  us  now  consider  the  case  (ii)  where  /is  an  internal  face  of  A  with  i  islands,  so 
that  /has  two  sides  that  are  not  directly  connected  to  each  other.  More  precisely,  each 
side  of /is  a  topological  disk  with  i  open  disks  removed,  so 

Xif)  =  2[x(D2)  -  ,-x(D2)]  =  2(1-/). 

Once  again,  to  obtain  the  boundary  dC'  of  the  new  cell(s)  we  have  to  delete  i  +  l 
cylinders  from  dC  (one  for  each  component  of  the  boundary  of/)  and  glue  it  to  /along 
the  2(j-l-l)  "seam"  circles.  The  Euler  characteristic  of  dC  with  /  + 1  cylinders  removed 
is  xi^C).    "Gluing"  it  to  /along  the  2(i  +  1)  circles,  we  obtain 

XOC)  =  x(3C)  +  x(/)  =  x(3C)  +  2(l-/). 

Recalling  the  expression  for  the  Euler  characteristic  of  the  components  of  dC  and  dC 
met  by  /,  we  obtain  2(c'  -  g')  =  2(c  -  g)  +  2(1  -  i),  so  that  8^  =  8c  +  /  -  1.  Here 
again,  c'  is  either  2  or  1  depending  on  whether  /splits  C  into  two  cells  or  not.  (Clearly, 
if  C  is  split  in  two,  c'  =  2.  Otherwise,  consider  the  connected  component  K  of  dC  con- 
taining one  side  of/.  If  it  includes  the  opposite  side  of/,  c'  =  l,  as  asserted.  If  it  does 
not,  observe  that  AT  is  a  compact  connected  oriented  two-manifold  (without  boundary); 
thus  its  removal  disconnects  R''  into  an  interior  and  an  exterior  components  (see,  for 
example,  [Mo,  Chap.  26]).  As  K  does  not  include  the  other  side  of/,  two  points  lying 
sufficiently  close  to  each  other  on  opposite  sides  of /are  separated  by  KCdC ,  making  it 
impossible  for  both  of  them  to  lie  in  the  same  cell,  and  thus  contradicting  the  assump- 
tion that  C  was  not  split  in  two.)  On  the  other  hand,  c  can  be  as  high  as  j  +  1  or  as  low 
as  I  (/has  one  outer  boundary  and  /  islands;  the  extreme  cases  correspond  to  all  of  them 
lying  in  different  components  of  dC  or  all  lying  in  the  same  component).  Hence 
—  J  ^  8c  S  1  and  —  1  s  8g  ^  i.  Moreover,  if  /does  not  split  C  into  two  cells,  c'  =  l 
and  the  above  relation  becomes  8^  =  j  — c.  In  particular,  we  have  shown  that  each 
cutting-but-not-splitting  internal  face  /must  be  of  one  of  the  following  two  types: 

(1)  /  =  0  and  8^  =  — 1,  or 

(2)  i>0. 
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To  summarize  the  cases:  introduction  of  a  face  /  (in  Q  (A))  either  reduces  the  total 
genus  by  one  or  increases  it  by  (at  most)  the  number  of  its  islands  (if  /  is  a  boundary 
face  of  Q(A),  it  may  also  cause  an  additional  increase  in  the  total  genus  by  as  much  as 
the  number  of  components  of  doutf-  3A  less  1).  However,  the  number  of  islands  over 
all  faces  of  Q(A)  is  easily  seen  to  be  bounded  by  the  total  number  f^  of  segments  in 
Q(A)  (since  a  segment  cannot  appear  in  two  islands)  and  the  number  of  components  of 
5oui/~3A  over  all  faces /in  Q(A)  is  bounded  by  Ir^-   Hence: 

Theorem  A.l:  The  change  in  the  total  genus  of  all  cell  boundaries  in  an  arrangement  A 
of  triangles,  caused  by  adding  a  triangle  A  to  A,  is  between  3f ^  and  —k,  where  k  is  the 
number  of  internal  faces  of  Q(A),  having  on  islands,  whose  introduction  does  not  split  a 
cell  of  A  into  two  subcells,  and  t^  is  the  number  of  segments  in  Q(A). 

Thus,  summing  these  changes  in  an  incremental  construction  of  AiG),  as  in  the 
proof  of  the  Combination  Lemma  in  Section  2,  we  obtain: 

Proposition  A. 2:  In  an  incremental  construction  of  an  arrangement  of  triangles,  as  car- 
ried out  in  the  proof  of  the  Combination  Lemma,  the  total  number  of  faces  that  cut  cells 
without  splitting  them  is  at  most  0(1),  where  t  is  the  number  of  pairs  of  triangles  that 
have  non-empty  intersection. 

Proof:  We  will  count  the  number  of  such  faces  in  a  construction  that  starts  with  no  trian- 
gles and  proceeds  to  add  first  all  red  and  then  all  blue  triangles,  one  at  a  time.  This 
clearly  provides  an  overestimate  on  the  desired  quantity.  The  total  genus  of  the  boun- 
daries of  all  cells  is  always  non-negative  and  the  genus  of  an  empty  set  (the  boundary  of 
R^— the  only  cell  in  an  arrangement  of  no  triangles)  is  zero.  It  is  decreased  by  one  by 
the  introduction  of  each  type  (1)  internal  face.  On  the  other  hand,  the  previous  theorem 
implies  that  the  overall  increase  in  the  total  genus  is  proportional  to  the  total  number  of 
segments  in  all  arrangements  Q(A),  which  is  to  say,  0(t).  Hence  the  total  number  of 
faces  for  which  8^  =  — 1  is  0(t).  As  islands  cannot  be  shared  among  faces  of  Q(A), 
there  are  at  most  t^  faces  of  type  (2)  in  Q(A),  for  a  grand  total  of  0(t).  The  conclu- 
sion follows.   □ 

Theorem  A.3:  The  genus  summed  over  the  boundaries  of  all  cells  in  an  arrangement  of 
triangles  in  which  t  pairs  intersect  is  at  most  0(0- 

Proof:  Immediate  from  Theorem  A.l. 

A.2  Euler's  relationship  for  cell  boundaries 

In  this  section  we  will  justify  our  practice  of  bounding  the  complexity  of  a  cell  by 
the  number  of  its  faces.  Consider  a  single  connected  component  AT  of  a  cell  boundary. 
By  using  the  puffy  model,  we  ensure  that  AT  is  a  compact  connected  orientable  2- 
manifold,  and  thus  we  can  apply  standard  topological  techniques  to  relate  the  number  of 
faces,  edges,  and  vertices  in  K  to  its  Euler  characteristic  and  genus  (see,  for  example, 
[GP]).  Since  faces  of  K  are  not  necessarily  simply  connected,  we  must  account  for  this 
as  well.  Let  vg,  e^,  fK  be  the  number  of  vertices,  edges,  and  faces  of  K,  respectively. 
If  all  faces  of  K  were  simply  connected,  the  definition  of  the  Euler  characteristic  x(A^)  of 
K  would  yield 

VR-eK-^fK  =  X(K). 

Consider  a  non-simply  connected  face— its  boundary  consists  of  a  single  outer  com- 
ponent and  zero  or  more  islands.    Notice  that  introduction  of  /  cuts  transforms  an  /- 
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island  face  into  a  simply-connected  one,  where  a  cut  is  a  simple  arc  connecting  a  vertex 
of  an  island  to  a  vertex  of  the  outer  boundary  of  the  face.  Thus  all  non-simply- 
connected  faces  of  K  can  be  eliminated  by  increasing  the  number  of  edges  in  K  by  the 
total  number  ijc  of  'islands'  on  the  faces  of  K.  Hence 

VK  -  (eic  +  Ik)  +fK  =  X(K)    or    vk  +  fg  =  x(^)  +  ^a:  +  'a:.  (**) 

Notice  that  all  vertices  of  K  besides  triangle  corners  have  degree  three  or  more,  while 
triangle  corners  have  degree  two.  Denoting  the  number  of  such  corners  on  K  by  c^  and 
summing  the  degree  over  all  vertices,  we  obtain 

2ck  +  3(va:  -  cat)  ^  lex    or     Zvk  ^  le^  +  ck-  (***) 

Conditions  (**)  and  (***)  together  imply 

eK^CK  +  2fK  -  3x(K)     and     vic^Cfc  +  Vk  -  2x(^) • 

Recall  that  the  Euler  characteristic  of  a  compact  connected  orientable  manifold  is 
2(1  —  gK)y  where  g)^  is  the  genus  of  K,  yielding 

^a:  <  CAT  +  3/jsr  +  dgK, 

and  a  similar  bound  for  vj^ .  Summing  these  inequalities  over  all  connected  components 
K  of  dC  (with  C  being  a  fixed  cell  or  any  collection  of  cells),  we  obtain 

«dc  -  cac  +  3/ac  +  6^ac  =  ©(/ac  +  n  +  O 
and 

Vdc  ^  cac  +  2/ac  +  4gac  =  C)(/ac+"+0, 

where  g^c  and  c^c  refer  to  the  total  genus  of  all  components  of  bC  (in  the  sense  of  the 
last  section)  and  the  total  number  of  triangle  corners  on  these  components,  respectively. 
The  last  equality  follows  from  Theorem  A. 3  and  the  observation  that  there  are  3n  trian- 
gle corners  in  the  whole  arrangement.   Thus  we  have  shown: 

Theorem  A. 4:  The  total  combinatorial  complexity  of  a  collection  of  cells  in  an  arrange- 
ment of  n  triangles  in  space  is  0(F  +  r  +  n),  where  F  is  the  total  number  of  faces  in  the 
boundaries  of  these  cells  and  t  is  the  number  of  intersecting  pairs  of  triangles.  Q 

Appendix  B.  Other  combination  lemmas— an  example. 

In  this  appendix  we  demonstrate  the  strength  of  the  technique  that  we  used  in  prov- 
ing the  Combination  Lemma  (Lemma  2.1)  by  employing  it  to  obtain  a  simple  proof  of  a 
variant  of  the  combination  lemma  of  [EGS2]  for  faces  in  arrangements  of  segments  in 

Lemma  B.l:  Consider  two  arrangements  of  segments,  red  and  blue,  with  a  total  of  n 
segments.     Assume    that    no    two    segments    intersect    in    more    than    a    point.     Let 

B=  {B\,  .  .  .  ,5j}  (resp.    11=  {R\ /?,})  be  a  family  of  faces  in  the  red  (resp.  blue) 

arrangement.    Let  P  be  a  set  of  k  points  {p\ pt},  such  that  each  point  p,  lies  in 

fiji  ^  ^f(  for  some  unique  5,-  and  f,.  For  each  /  =  1,  .  .  .  ,;t,  let  E,  denote  the  face  of  the 
combined  arrangement  (i.e.,  the  connected  component  of  Bj.DR,.)  containing  /?,.  Then 
the  complexity  of  the  "purple"  family  {£,}  is  at  most  ^  +  p  +  0{k+n),  where  3,  p  are 
the  complexities  of  B,  V.,  respectively,  and  the  complexity  of  a  face  containing  two  or 
more  points  of  P  is  still  counted  only  once. 


Triangles  in  space  -  34  -  October  21,  1988 


Proof.  Let  tir  and  tif,  be  the  number  of  red  and  blue  segments,  respectively.  We  follow 
the  general  approach  of  the  proof  of  Lemma  2.1,  but  the  analysis  here  is  considerably 
simpler.  That  is,  starting  with  the  red  arrangement,  we  will  incrementally  superimpose 
the  blue  segments  on  it,  thereby  trimming  and  subdividing  the  red  faces  in  11  until  they 
assume  their  final  "purple"  shape.  In  the  process  we  will  bound  the  number  of  addi- 
tional red  edges  that  are  created  on  purple  boundaries.  Repeating  the  same  process, 
starting  with  the  blue  family  B  and  incrementally  transforming  it  to  the  purple  family, 
will  provide  an  estimate  on  the  number  of  additional  blue  edges  bounding  the  E,'s. 
Together,  these  bounds  give  an  upper  bound  on  how  much  the  complexity  of  {E,} 
exceeds  p  +  p. 

We  thus  begin  with  the  red  family  TZ  and  add  blue  segments,  one  by  one.  At  any 
step  during  this  process,  we  refer  to  the  polygonal  regions  of  the  current  planar  map 
containing  points  of  P  as  'currently  purple'.  Consider  the  next  blue  segment  e  to  be 
added  to  the  arrangement.  We  will  refer  to  a  maximal  segment  of  the  intersection  of  e 
with  a  currently  purple  face  as  a  blue  fragment.  Clearly  the  next  currently  purple  family 
may  be  obtained  by  adding  all  blue  fragments  of  e  to  the  currently  purple  regions  and 
then  considering  the  regions  of  the  resulting  planar  map  which  contain  points  of  P. 
Consider  adding  blue  fragments  to  the  current  arrangement  one  at  a  time.  Let  fe  be  a 
blue  fragment  (contained  in  some  currently  purple  face  Q)  being  added  in  the  current 
step.  By  construction,  the  relative  interior  of  b  does  not  meet  any  currently  purple  face 
boundaries.   Therefore,  b  can  be  classified  into  exactly  one  of  the  following  five  classes: 

(i)  b  cuts  Q  into  two  subfaces,  each  containing  points  of  P; 

(ii)  b  cuts  off  a  portion  of  Q  containing  no  point  of  P\ 

(iii)  b  connects  two  components  of  3Q,  without  fully  cutting  Q; 

(iv)  b  meets  dQ  at  only  one  of  its  endpoints; 

(v)  Z?  is  a  full  blue  segment  lying  entirely  in  the  interior  of  Q. 

(Note  that  the  class  of  a  fragment  is  determined  at  the  time  of  insertion.)  We  now 
proceed  to  estimate  the  total  number  of  blue  fragments  in  each  class  over  the  course  of 
construction  and  the  amount  of  additional  red  complexity  that  they  create.  In  case  (v), 
introduction  of  b  clearly  has  no  effect  on  the  red  complexity  of  currently  purple  faces, 
so  let  us  restrict  our  attention  to  blue  fragments  of  types  (i)-(iv).  Each  fragment  of  type 
(i)  or  (iii)  creates  at  most  two  new  red  edges  (as  it  meets  dQ  exactly  twice).  There  are 
at  most  k-\  type  (i)  fragments,  as  each  such  fragment  further  refines  the  partition  of  P 
into  sets  of  points  lying  in  distinct  currently  purple  cells.  As  each  type  (iii)  fragment 
connects  two  'islands'  of  a  currently  purple  cell  and  no  segment  (either  blue  or  red)  can 
lie  in  two  islands,  type  (iii)  fragments  occur  at  most  n-\  times.  In  particular,  there  are 
at  most  /If— 1  type  (iii)  fragments,  connecting  two  points  of  red  segments  and  thus 
increasing  the  red  complexity  by  2,  and  at  most  n^  fragments  of  type  (iii)  linking  a  point 
on  a  previously  inserted  blue  segment  to  one  on  a  red  segment,  each  increasing  the  red 
complexity  by  1.  The  fragments  of  type  (iii)  that  connect  two  points  of  previously 
inserted  blue  segments  do  not  affect  the  red  complexity.  As  for  type  (ii)  fragments, 
note  that  the  endpoints  of  b  lie  on  different  edges  of  dQ,  as  6  is  a  line  segment  and  no 
two  segments  overlap.  Each  such  edge  is  cut  by  the  endpoint  of  b  into  two  subedges, 
one  of  which  lies  outside  of  the  newly  truncated  currently  purple  cell,  and  thus  no  addi- 
tional red  edges  are  created.  Turning  to  type  (iv)  fragments,  each  of  them  creates  at 
most  two  new  red  edges  out  of  the  edge  that  it  meets,  thereby  increasing  the  red 
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complexity  by  no  more  than  one.  However,  the  endpoint  of  b  that  does  not  lie  on  BQ 
must  be,  by  construction,  an  endpoint  of  a  blue  segment,  so  that  there  are  at  most  2nb 
type  (iv)  blue  fragments.  Thus,  once  all  blue  segments  have  been  added,  the  overall 
increase  in  the  number  of  red  edges  in  all  purple  regions  is  at  most 

2(k-l)  +  2(nr-\)  +  ni,  +  2nb  =  2A:  +  2n,  + 3/1^-4. 

Repeating  the  analysis  for  the  increase  in  blue  complexity,  we  conclude  that  the  total 
complexity  of  all  purple  faces  is  at  most 

P  +  p  +  4A:  +  5rt-8, 

as  asserted.  □ 
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